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THE STABILITY OF THE LAMINAR BOUNDARY LAYER 
IN A COMPRESSIBLE FLUID 


By Lester Lees 
SUMMABY 


The present paper is a continuation of a theoretical invest! - 
^tion of the stability of the laminar boundary layei’ in a com- 
pressible fluid. An approximate estimate for the critical 

Reynolds number Eg , or stability limit, is obtained in terms 

of the distribution of the Islnomatic viscosity and the product of 


the mean density p* and mean vorticity aci'oss the boundary 

layer. With the help of this estimate for Rg , it is shcrwn 
■ *^^rjin 

that vithdraving heat from the flxxid through the solid surface 
increases stabilizes the flov, as compared with the 

flow over an insulated surface at the same 'Mach number. Conduction 
of heat to the fluid through the solid surface has exactly the 
opposite effect. The value of R© for the insulated surface 

•^^min 

decreases as the Mach number increases for the case of a uniforia 
velocity . These general conclixsions are supplemented 
y detailed calculations of the curves of wave number (inverse 
^ve length) against Reynolds number for the neutral disturbances 
for 10 representative cases of insulated and noninsulatod surfaces. 


So far as laminar stability is concerned, an i jTip ortant dif- 
ference e^sts between the case of a subsonic and. supersonic froe- 
stroam velocity outside the boundary layer. The neutral boundary- 
layer disturbances that are significant for laminar stability die 
out oponontially with distance from the solid surface; therefore 
ho phase velocity c* of tho^ disturbances is subsonic relative 


* - 


c* < a^*. 


•whore 


a * 
o 


■fco the freo-streaui *7olocity lu* - or u 
is the local sonic velocity. When ^ < 1 (;diere is 

free-stream Mach number), it follows that 0 < < c* ; and any 
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laminar toundary-layer flov is iiltimately unstalDle at sufficiently 
hi{^ Eeynolds numbers because of the destabilizing action of vis- 
cosity near the solid surface^ as e3(plained by Prandtl for the 


incompressible fluid, VUien Mq > 1 ^ 

If the quantity |-^ |p* )j 
Ldy* \ 

the rate at -which energy passes from 


ho-wever, ^ >1 - > 0. 

is large enou,^ negative.ly, 
the disturbance to the mean 


flov, which is proportional to 



can 


al-ways bo large enou^i to counterbalance "the rate at -vdiich energy 
passes from the moan flov to tlae distm’bance because of the desta- 
bilizing action of viscosity near, the solid sur'face. In that case 
only damped disturbances exist and the laminar boundary layer is 
completely stable at all Rejnolds numbers. This condition occurs 
when the rate at \diich heat is tvi-bhdra-wn from the fluid throu^ 
the solid surface reaches or exceeds a critical value that depends 
only on the Mach numbez- and the propei-ties of the gas. Calcula- 
tions show that for Mq > 3 (approx.) the laminar boundary -layer 


flow for thermal equilibrium - viiere the heat conduction tlzrougja 
the solid sui'face balances the heat radiated from the surface - is 
completely stable at all Eeyno.lds numbers under free-fliglit conditions 
if the free-stream velocity is unifomi. 


The results of the analysis of the stability of the laminar 
boundary layor must bo applied -with care to discussions of transi- 
tion; however, \d.thdrs-Ming heat from the fluid throu^ the solid 
surface, for example, not only increases Eg but also 

'^^min 

decreases the initial rate of aJ3g)lifi cation of the aelf-oxcitod 


disturbances, which is roughly proportional to 


cr 


Thus, 


nun 


cr. 


tho effect of the thenaal conditions at the solid surface on tlie 
transition Eeynolds number ^ 9 ^^ -s similar to the effect on Eg 

A comparison between -this conclusion and experimentai. investigations 
of the effect of surface heating on transition at low speeds shows 
that the results of the present paper give the proper dj.rection of 
this effect. 


nun 


The extension of tho results of the stability analysis to 
laminar boundary-layer gas flows -vdth a pressiire gradient in tho 
direction of the free s-tream is discussod. 
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INTROD'UGTION 


By blie theoretical studies of Heisenberg, Tollmien, Schlichting, 
and Lin (references 1 to 5) careful esroeriiiiental investi- 

gations of Liepiaann (reference 6) and H. L. Drydon and liis asso- 
ciates (reference 7), ib has been definitely established tliat the 
flow in the laaiinar boundary layer of a viscous hesmogeneous incom- 
pressible flvdd is unstable above a certain characteristic critical 
Reynolds nuiaber. Wlien the level of the disturbances in the free 
stream is low, as in most cases of teclin5.cal interest, this inherent 
instability of the laminar motion at sufficiently hi^ Reynolds 
numbers is responsible for the ultimate transition to turbulent 
flow in the boundary layer. The steady laminar boundary-layer flow 
would always represent a possible solrition of the steady eciuations 
of motion, but this steady flow is in a state of unstable dynamic 
eq.uilibrium above the critical Reyno3.ds number. Self -excited dis- 
turbances (Tollmlen waves) appear in the flow, and these disturb- 
ances grow large onou^ eventually to destroy the laminar motion. 

The question naturally arises as to how the phenomena of 
laminar instability and transition to turbulent flow are modified 
when the fluid, velocities and temperature variations in the boimdary 
layer are large enoui^ so that the coBipressibility and conductivity 
of the fluid can no longer be neglected. The present paper repre- 
sents the second phase of a theoretical investigation of the sta- 
bility of the laminar boundary-layer flow of a gas, in vMch the 
compressibility and heat conductivity of the gas as v’ell as its 
viscosity, are taken into account. The first part of this work 
was presented in reference 8. The objects of this investigation 
are (l) to determine liow the stabii.ity of the lamincar boundary 
layer is affected by the free -stream Mach number and the thermal 
conditions at the solid boundary and (2) to obtain a better under- 
standing of the piprsical basis for the instability of laminar gas 
flows. In this sense, the present study is an extension of the 
Tollmien-Schlichtlng analysis of the stability of the laminar flow 
of an incompressible fluid, but the investigation is also concerned 
with the general question of bo\mdary-layor disturbances in a 
con5>roBsiblo fluid and tlieir possible interactions -s-dth the main 
external flow. 


SnffiOLS I 

I 

I 

With minor exceptions the symbols used in tills paper are the I 

same as those introduced in reference 8. Phjrsical quantities are | 
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denoted "by an asterisk, or star, •vdaereaa the corresponding non- 
diaensional quantities are unstarred. A "bar over a quantity denotes 
mean value; a prime denotes a fluctuation; the suhscript o denotes 
free -stream values at the "edge*' of the boundary layer; the sub- 
script 1 denotes values at the solid surface; and the sub- 
script c denotes values at the inner “critical layer”, T-iiore 

the phase velocity of the disturbance equals the mean flow velocity. 
The free-stream values are the characteristic measures for all non- 
dimensional quantities. The charactez-istic length measure is the 
boiUidary-layer thickness 6, except •vdioi'o otheia-dse indicated, 
llote that in order to confoim vith standai’d notation, the symbol & 
for bcimdary-layer thickness is unstazrred., vdiereas the syonbols 5* 
and 0 are zised for boundaz'y-layer displacement thickness and 
boundary-].ayer momentum thickness, respectively. 

X* distance along surface 

distance normal to surface 

t* time 


u* 



V* 


V*’ 

iD = 

V 


P* 

P* 

T-x- 

T* 


component of velocity in -direction 


coDTOonent of velocity in y^-d.iroct:lon 


stream function for mean flow 
density of gas 
pressure of gas 
temperature of gas 
laminar sheaz* stress 

ordinary coefficient of viscosity of gas 
kinematic viscosity of gas 
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P 

E«- 

7 

c* 

x->;- 

5 

5* 


thermal conductivity of gas 
specific heat at constant volume 

specific heat at constant pressure 
gas constant per gram 

ratio of specific heats ^'^p/^v)^ 1.^5 ^or air 

complex phase velocity of boundary-layer disturbance 
wave length of boundary-layer distirrbance 
boundary-layer thi cameos 

/r 

boundary-layer displacement tliicloaess [I (l - pw)dy* 


boundary-layer momentum thickness 


pv(l - w)dy^ 


a* 


a - 


Ctg = 


E 


23T 

X*/5 

2n 

X*/G 


1 -rave number of boundary-layer disturbance (2n/X*) 


Eeynolds number 


"po* Uo'» 6\ 




/ 


= 


Pq* Uq* 0 






Mach nvmibor 


V 


\/7E^ Tq* 
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C 


Prandtl number jcp 


) 


1 . PRSL.BIIWAEY COWCIDERATIOWS 


In the first phase of this investigation (reference 8 ) the 
stability of the laminar boundary-layer flov of a gas is analyzed 
by the method of small perturbations, which I'jas already so suc- 
cessfully utilized for the study of the stability of the lamlnaz' 
flow of an incompressible fluid. (See reference 5 .) by this 
method a nonsteady gfis flow is investigated in \diich all physical 
q.unntities differ from their values in a given steady gas flov 
by small perturbations that are functions of the time and the space 
coordinates. This nonsteady flow miJist satisfy the complete gas- 
dynamic equations of motion and the same boundary conditions as 
the given steady flow. The question is whether the nonsteady flow 
damps to the steady flow, osci.llates about it, or diverges from it 
■vrit?! time - that is, 'sdiether the small perturbations are dairped, 
neutral, or self-excited disturbances in time, and thvis whether 
the given stoad;;,’' gas flow is stable or unstable. The analysis is 
particularly concerned with the conditions for the existence of 
neutral disturbances, idiich mark the transition fz’om stable to 
unstable flow and define the minimum critical Reynolds nimiber. 

In order to bring out some of the principal features of the 
stability problem ■'.ri.thozit becoming involved in hopeless mathe- 
matical cesmpli cations, the solid boundary is taken as two dimen- 
sional and of negligible curvature and the boundary-layer fd.ow is 
regarded as plane and essentially parallel; that is, the velocity 
component in the direction noimal to the surface is negligible and 
the velocity coiaponent parallel to the sxirface is a function mainly 
of the distance normal to the surface. The small disturbances, 
which are also two diDiensional, are analyzed into Fourier com- 
ponents, or noriiial modes, periodic in the direction of the free 
stream; and the aaiplitude of each one of these partial oscillations 
is a function of the distance normal to the solid surface, that 


In the study of the stability of the laminax’ boundary layer, 
it will be seen that only the local properties of the "parallel” 
flow are significant. To include the variation of the mean velocity 
in the direction of the free stream or tho velocity component noarmal 


is, 
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to the solid boundary in the problem v’Duld lead only to hij^er order 
terms in the differential equations governing the disturbajicea, 
since both of these x'actors are invorseljr proportional to the local 
Reynolds number based on the bovindary-layer thickness. (See,, for 
exar!q)lo, reference 2.) By a careful analysis, Pretsch has shown 
that even with a pressure gradient in the direction of the free 
stream the local mean-velocity distribution alone determines the 
stability characteristics of the local boundary-layer flow at 
large Reynolds numbers (reference 9 ) • Such a statement applies 
only to the stability of the f.low vdthin the borada.ry layer. For 
the interaction between the boundary ].ayer and a main “external" 
supersonic flow, for example, it is obviously the variation in 
bouTidary-layer thickness and mean velocity along the surface that 
is significant. (See reference 10.) 

The aforementioned considerations also lead quite naturally 
to the stud^ of individual partial oscillations of the 

form f(y) for which the differential equations of 

distmbanco do not contain x and t explicitly , These partial 
oscillations ai-e ideally suited for the study of instability, for 
in ordei- to show that a flow is unstable it is unnecessary to 
consider the most gone-ral poss3.blo disturbance; in fact, the 
simplest will suffice. It is only necessary to show that a 
particular disturbance satisfying the squatio.ns of motion and the 
boiuidary conditions is self-excited or, in this case, that the 
imaginary part of the comp.lex phase ve3.ocity c is positive. 

In reference 8 the differential equations governing one 
normal mode of the disturbances in the 3nminar boundary layer of 
a gas were derived and studied very thorou{dil 5 ’". The conq)lete set 
of solutions of the distm’bance equations ’.ras obtained and the 
physical boundary conditions that these solutions satisfy were 
investigated. It was found that the final relation between the 
values of c, a, and R that determines the possible neutral 
disturbances (limj.ts of stability) is of the same form in the 
compressible fluid as in the incon^re'ssible fluid, to a first 
approximation. The basis for this result is the fact that for 
Royno.lds numbers of ' the order of those encountered in most aero- 
djmamic problems, the temperature disturbances have only a negligible 
effect on those parti cul.ar vc.locity solutions of the disturbance 
equations that depend primarily on the viscosity (viscous solu- 
tions).- To a first approximation, those viscous solutions there- 
fore do not depend directly on the heat conductivity and are of 
the same form as in the incompressible flviid, except that they 
involve the Rojaiolds number based on the kinematic viscosity near 
the solid boundary (where the viscous forces are important) rather 
than in the free stream. In this first approxiiHation, the second 
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viscosity coofficiont^ viiich. is a measure of the dependence of the 
pressure on the rate of change of density^ does not affoct the sta- 
hility of the laminar hoimdary layer, Fi'om these results it vas 
inferred that at large Roynolds numbers the iniluence of the viscous 
forces on the stability is essentially the same as in an incom- 
pressible fluid. This infei^ance is borne out by the results of the 
present paper. 


The influence of the inertial forces on the stability of the 
laminar boundary layer is reflected in the behavior of the asymp- 
totic inviscid solutions of the disturbance eguations^ idiich are 
independent of Reynolds number in first approx 3 .mation. The results 
obtained in reforenco 8 show that the behavior of the inertial 
forces is dominated by the distribution of the product of the moan 


dw 

densitj' and moan vorticlty p — across tho boundary layer. (Tho 

^ d / dv\ 

gradient of this quantity, or — (p — ), -vdiich p 3 .ays the aamc role 

dy \ dy / 

as the gradient of the vorticity in the case of an incomprossiblo 
fl’uid, is a measure of the rato at "which the :>:-momentum of the 
thin layer of fluid near the critical layer (tdiero v =* c) 
increases, or decroasos, becavtso of the transport of momentum by 
the disturbance.) In order to clarify the bohavior of the inortial 
forces, tho limiting case of an inviscid fl.uld (E— is ctiidiod 
in detail in rofei’enco 8. Tho folloTring goneral cr iter ions are 

obtained: (l) If tho quantitjr l(p^\ vanishes for somo value 

of w > 1 - i-. then neutral and self -oxci tod subsonic disturb- 
ancoa exist and the inviscid compressible flov Is uristablo. 


(2) If tho quantity 


/cb.’ 
■ P 


3 ^ 


does not vanish for some valuo 


01 




thon all subsonic disturbances of finite tovo 


length aro danced and the inviscid coitprossiblo flov is stable. 
(Outside the boundary layer, tho relative velocity betveon the moan 
flow and tho x-component of tho phase velocity of a subsonic dis- 
turbance is loss than the mean sonic velocity , Tho magnitude of 
such a distrcrbance dios out o>5)onentially with distance from the 
solid suT’facG.) ( 3 ) In goneral, a dioturbanco gains energy from 
0 /• dv\ 

the mean flow if ~(p— - | is positive at the critical layer 

<iy\<3.37 

(where - c) and loses energy to tho moan flow if 
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The general stability criterions for inviscid compressible 
flov give some insi^it into the effect of the inertial forces on 
the stability, but they cannot bo taJcen over bodily to the real 
compressible fluid. Of course, if a flow is unstable In the . 
Ijjuiting case of an infinite Reynolds number, the flow is unstable 
for a cortain finite range of Reynolds number. A coii5)resaible flow 
that is stable when R — co, however, is not necessarily stable 
at all finite Reynolds numbers vixon the effect of viscosity is 
taken into account. One of the objects of the present paper is 
to settle this question. 


On the basis of the stability criterions obtained in refer- 
ence 8, some general statements were made concerning the effect of 
thermal conditions at the solid boundary on the stability of laminar 
boiuidary-layer flow. It is concluded, from physical reasoning and 


a study of the equations of mean motion tliat the 


vanishes for some value of w > 0 



quantity 



that is, if 


heat is a dded to the fluid tliroug^i the solid surface or if the 
suri-acQ is insulated. If > 0 and is sufficiently large, 

that is, if heat is Xidthdrawn from the fludd throu^ the solid 

surface at a sufficient rate, tlie quantity ] never vanishes, 

dv dv f 


Thus , when 


(-) 

Wi- ■ 


dy dy / 

the laminar boundary-layer flow is desta- 


bilized by the action of the inertial forces but stabilized 
throxigh the increase of kinematic viscosity near the solid surface. . 

Vhen ■bbe reverse is true. The question of -vdiich of 


these effects is predominant can be answered only by further study 
of the stability problem in a real congsressible fluid. 


In the present paper this investigation is continued along the 
f olio v/ing, lines; . ■ • 

(1) A study is made of how the general criterions for insta- 
bility in an inviscid cctapressible fluid are modified by the 
introduction of a small viscosity (stability at very laz'ge 
Reynolds numbers) . • • 
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(2) The conditions for the existence of neutral disturb- 
ance s at large Reynolds nuoiber are examined (study of asymp- 
totic form of relation between eigen-values of c, a, and R). 


(3) A relatively simple expression for the approximate 
value of the minimum critical Reynolds number is derived-, 
this expression involves the local distribution of mean 
velocity and mean temperature across the boundary layer. This 
approximation will serve as a criterion from which the effect 
of the free- stream Mach nxmber and theimial conditions at the 
solid surface on the stability of laminar boundary- layer flow 
is readily evaluated. The question of the relative influence 


of the kinematic viscosity and the distribution of 
stability would then be settled. 



on 


(h) Tho energy balance for smeill disturbances in the roal 
compressible fluid is considered in an attempt to clarify the 
physical basis for the instability of laminar gas flows. 

(5) In order to supplement the investigations outlined 
in the four preceding paragraphs, detailed calcixlations are 
made of the limits of stability, or the curve of a against R 
for the neutral disturbances for several representative cases 
of insulated and noninsulated surfaces. The results of the 
calculations are presented in figures 1 to 8 and tables I 
to IV . The method of compxxtation of the stability limits is 
briefly outlined in reference 8 , although the calculations 
were not carried out in that paper. 

In the present investigation the work of Heisenberg (refer- 
ence 1) and Lin (reference 5) on the stability of a real incom- 
pressible fluid is naturally an indispensable guide. In fact, the 
methods utilized in the present study are analogous to those 
developed for an incompressible fluid. 

The present paper is concerned only with the subsonic disturb- 
ances. The amplitude of the subsonic disturbance dies out rapidly 
with distance from the solid boundary. In other words, the neutral 
subsonic disturbance is an "eigen-osclllatlon" confined mainly to 
the boundary layer and exists only for discrete eigen-values of c, 
a, and R that determine the limits of stability of laminar 
boundary-layer flow. Disturbances classified in reference 8 as 
neutral ’’supersonic," that is, dlstvirbances such that the relative 
velocity between the x-component of the phase velocity of such a 
disturbance and the free- stream velocity is greater than the local 
mean sound speed in the free stream, are actually progressive sound 
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■waves thai: 5japlnge obliquely on the boundary layer and are reflected 
with change of amplitude. For disturbances of tlvia type •fclie -wave 
length and phase velocity are obviously completely arbitrary (eigen- 
values are continuous), and these disturbances have no sigaificance 
for boundary-layer stability, 

Men the free-stream velocity is supersonic > 1^, the 

subeon3.c boundary-layer disturbances must satisfy the requirement 

that Uq* - c* < ^ or c > 1 - i- (for < 1, c > O) . Now, 

^ 0 

by analoQT ■vdth the case of an incon 5 >resBible fluid it is to be 
espected that for values of c greater than some critical value of c^, 

say, all subsonic disturbances are daiaped. 'Thus, when Mq > 1, 

there is the possibility that for certain mean velocity -temperature 
distx-ibutions across the boxjndary layer, neutral or self-excited 
distva'bances satisfying the diffei-entiaD. equations of motion, the 

boundary conditions, and, also, the physical requirement tliat c > 1 - 

^o 

cannot be found. In that event, the laminar boundary flo'vr Is stable 
at all Eeynolds numbers. This Interesting possibility is investi- 
gated in the present paper. 


2. CALCULATTON OF TEE LBOTS OF STABILin OF THE LAlyUK'AR 
BOUNDARY LAYER EJ A VISCOUS COEDUCTIVS GAS 

order that tho coDplete system of solutions of the differ- 
ential equations for tho propagation of small disturbances in the 
laminar boundary layer shfill satisfy the physical boundary condi- 
tions, tho phase velocity must depond on the wave length, the 
Reynolds nvimber, and the Mach number in a manner that is determined 
entii’oljr by tho local distribution of mean velocity and mean tempera- 
ture across the Doundary layer. In other words, tho only possible 
suDGonic disturbances in tho laminar boundary layer are those for 
Mich there exists a definite relation of the form (reference 8) 

c - c(a, R, (2.1) 

Since a, R, and are real quantities, the relation (2.1) is 

equivalent to the t\ro relations 
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(2.1a) 

°i = (2.1b) 

The curve R_. Mq^) = 0 ^or a = a^E, ^o^)) neutral 

disturbancea gives the limits of stability of the lamliaar boxmdary 
layer at a given value of the Maoh number. From this ciorve can be 
determined the va3,uo of the Reynolds numbei’ be3.ow vhich disturbances 
of all viave lengths are damped and above "/iiich self -excited disturb- 
ances of cei’tain i-ravo lengths appear in a given laminar boundary- 
layer flow. 

3h reference 6, it is shown that the relation (2.1) between 
the phase velocity and tlie wave 3.en.3th takes the following form; 


E(a, c, = F(z) 


( 2 . 2 ) 


In equation (2.2), F(z) is the Tietjens function (reference 11) 
defined by the relation 



and the quantity Hanlcel fimction of the first kind 

of order I/ 3 • The prime denotes differentiation with respect 
to y. The function c, x-hich depends only on the 
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asymptotic invlscicL solutions end <pg (section 4 of refer- 

ence 8) and not on the Reynolds nvouber, is defined as follows; 


(7l - 7c) ^o^) = 


<? n 9 ip * + P9 ip 

921 ‘+’ 22 * P9 pp 


Ti9ii‘ + Mo ^’^i’c9ii 

Ti - Mq2c2 


Ti92i’ •!- Mo^vri’ccpgi 


(2.5) 




9i2* + P9i2 


9p2* •^'^22 


where 


& = avjl - Mq2(i _ c)2 
i, J = 1, 2 


(2.6) 


and and are the coordinates of the solid surface and the 

”edge” of the houndai’y layer, respectively. 


The Tietjens function was carefully recalculated in reference 8, 

and the real and imaginary parts of the function $(z) — - 

1 - F(z) 

are plotted in figure 9* (The function 9(z) is found to he more 
suitable than F(z) for the act’jal calcu3.ation of the stability 
limits . ) 


The inviscid solutions and cpg were obtained as power 

series in or as follow® (section 6 of reference 8) : 
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OP 

cp^(y; c, = (v - c) ^ 


n=0 


9, 


00 

,(j'; c, Mq^) (v - c) ^ a^\2n+l(^"' 

rt=0 


vhere for n > 1 


h^Jy; 0 , = 


'iy ^ 


'J7^ 


T ', 2 
Mq 


.(v- c)^ 


ny 


ay 


uyi 


(tr - c)^ / o\ 

i^n-2(yj c> %‘-j^y 

I 

(2.9) 


and 


h = 1.0 
o 


and for n ^ 1 


>^.i(y: =,Mo®)= j 

J: 


py p 


T 2f 

— -V 


^yi 


(w- c)2 


iV ® 


and 


c, M^-) = 


r,7 r- 


T 


yy-, 


(w - c)2 


M. 




The lower limit in the integrcle is taken at the sijrface merely 
for convenience. Wien y > y^,, the path of integration must he 

taken below the point y = y_ in the complex y-plane. The power 

series in are then uniformly convergent for any finite value 

of a. 


Ilf 
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At the surface, the inviscid solutions ore readily evaluated 


’’-U “ 


”11' “ ’'1' 


-P21 - 0 


<p. 


21 




( 2 . 11 ) 


At tile "edge" of the boundary layer, the inviscid solutions 
most conveniently expressed as follows: 


are 


a%„(o, m/) 


n-u 


03 


m = (1 - 
^22 ^ 




n^O 


cp » 

^12 


(1 - c) 


1 - M 2(1 . c)2' 


(1 - c )2 


n=l 


( 2 . 12 ) 


^22’ - (1 - c) 


’1j^m£(^c)2‘ 

(1 - c )2 


-I CO 




n -0 


1? 
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where 




Ho = 1.0 




^an-lC*^^ ) ■■ 




1 - Mo^(l - c)' 

(1 -~ c )2 


1 - Mo^d - c)^ 
(1 - c)2 


n-1 


*'2n'(=^2> =< “oO 




(2.13) 


Ko = 1.0 


J 


With the aid of eq.\iations (2.11), tlie exjiression for ^(a, c, 
can he rewritten aa follo'\>?s: 


Wt_ • (922* •’' ^^22) 


E(a, c, Mo^) - ^ ^ 


'"l'(^ 22 ' c“ (^12’ ■'■ ^^12) 


(2.14) 


where 



c 


(2.15) 


The relation (2,2) between the phase velocity and the wave length 
is brouj^it into a forn more su5. table for the calcuJ.ation of the 
stability limits by making use of the fact that for real values 

of c the imaginary part of E (a, c, Mo^) is contributed largely 
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by the integral (The procedure to be folloved is 

identical vith that used by Lin In the lim?,ting case of 

the incoiapreB Bible fluid (reference purt IH) •) Define the 

function §(z) by the relation 



1 

1 - F(s) 


(2.16) 


Then, 





(1 H X) 


(u 4 iv) 
3. 1/ \{u^ 


where 


(2.17) 


u + 


^'^22 \ 

"^1 Vni;' 


(2..18) 


Equation (2.17) is equivalent to the two real relations 



(1 -y X)v 
(1 


(2.19) 


§3,(2) = (1 


X) 


u(l -V \u) Xv^ 

(1 + \u)- -i- X^'v^^ 


( 2 . 20 ) 


The real and imaginary parts of ^(s) are plotted against z in 
figivro 9. 


The dominant tenn in the imr.ginary part of the right-hand side 
of equation (2.l8), "idiich involves IC-, is exti'acted by means 

of straightforward algebraic transformations. Eolation (2.18) becomes 


17 


Vt ‘C 


u+ iv = 


T-, 






\]l - . c) 


TA a - a%c 


(1 - c)2 


2 X 00 




n=l 


rt=l 


2n+l 


a 




:n+l , 


( 2 . 21 ) 


Tdiore 




H 

CD 


and for n ^ 3 


and 


Tvlien c is real, 


l^n = %^-l - '"h 

(2.22a) 

\ = ^^n-2 - ^ 

(2.22b) 

’ c 


T sr I.P, K-, 

Ti 



for those values of a and c that occur in the stability calciilation^. (This approximation is 
Justified later in appendix A.) The imagii'tary part of the Integral 1^2^10, Mq^) is readily computed. 
It is found that 
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Now X(c) is cenerally quite Email, thgtrofore ^.,(z) caii be 

aken equal to v(c/ and <.52,(s) can be takcyn equal to u as a 
zerothiUpproxlmation. i^om equations (2.19) end (2.20), .dian c 



By equation (S.a!,) ,(o) ^ 


^ VV 

v'(l + X)3Vc y 


( 2 . 26 ) 


and the corresponding values of a 
(sliglitly transformed) by a method 


are obtained ;?rom equation (2.21) 
of successive approximations. 
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Thus, 


— (l-aHgj 


a = 


- 0)^ / 


1- <*%„)- A 


(1 - o)2 Vtlil '" '^V &1 “ "Sn+1_ 


(«.l) 




np2 




(2.27) 


vhere 


V, 'c 


"i - / ® i \ 

:r S.P.IKt ■; ) 

^J. \ ""I'V 


(The symbols and nov desJ.,^ate the real parts of the 

integrals Mjj. “and Nj^-. ) The Itei’ation process is begun by talcing 

a suitable initial value of , a on the ri^ht-liand. side of equa- 
tion (2.27). The- methods adopted for computing these intejyrals 
"When the mean velocity-temperature profile is Icnovai are described 
in appendixes A to C. 


For greater accuracy, the-, values of z and u for a given 
real value of c are computed by successive approximations. From 
equations (2.I9) and (2, 20), 




(1 t X)v 

(1 + X 


(2.28) 


Jn+1) ^ ^ (n+l)j^^(n4l)) 


(1 ^ aV ' 

Jl X) (1 


Xv 


-2 


(2.29) 


1 + Xu' 


,(n) 


The value of v ie alva7s approximated by relation (2.24). 

Curves of *v7ave number against Reynolds number for the neutral 
disturbance have been calculated for 10 representative cases 
(fig. 4), that is, insulated surface at Mach numbers of 0, O.50, 

0.70, 0.90, 1.10, and 1.30 and heat transfer across the solid surface 
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at a Mach mmher of O.7O with valueo of the ratio of surface teB5)era- 
ture to free-stream tenperatiire T^. O.80, O.90, am?. 1,25. 

(It is found more desii’ahle to base the nondlmensional \mre number and 
the Keynolds number on the momentum thickness G, "viiich is a dii'oct 
measure of the skin friction, rather tlian on tiie boundary-layer thick- 
ness 5 , which is somi^what indefinite.) , 

In figure 5 the jTiinimvim critical Ee;molds number , or the 

<=^min 

stability limit, is plotted against Mach number for the insulated 
surface; and in figure 6(a) Rq is plotted against T, for 

the cooled or heated surface at a liach number of 0 , 70 . The marked 
stabilizing infliuence of e. 'wlthdrauTul of heat from the flvd.d is 
clearly evident. Discussion of the physical significance of these 
numerical results is reserved until a.fter general criterions for 
the stability of the laminar bound-iry layer have been obtained. 


3 . DESTABILIZTOG ITIFIliRirCE OF VISCOSITY AT VERY LARGE EUTTOOLDS 
KUMEERS; EXTETISION OP HEISOI'BIRC- *S ORITSRION 


TO th :>3 comrresstble fluid 


The numerical calculation of the limits of stability for several 
particular cases gives some indication of the effects of free-stream 
Mach number and therm£‘l conditions at tho solid surface on the sta- 
bility of the laminar boundary .layer. It woirld be very desirable, 
however, to establish general criterions for Laminar instability. 

For tho incompressible fluid, Heisenberg has shown that tho influence 
of viscosity is generally destabilizing at very large Reynolds 
numbers (reference l) . His criterion can.be stated as follows: If 

a neutral disturbance of nonvanishing pliaso velocity and finite TOve 
length exists in an inviscid fluid (R-^oo) for a given moan velocity 
distribution, a disturbance of tho same wavo lonrctli is laistable, or 
self-excited, in the real fluid at very largo (but finite) Reynolds 
numbers . 


Tho same conclusion can bo drG•^^’n from Prandtl’s discussion of 
the enorgj’’ balance for small distimbanccs in tho lam3.nar boundary 
layer (reference 12) . 

Heisenberg’s criterion is established for subsonic disturbances 
in the laminar boundary layer of a compressible fluid by an argument 
quite similar to that which he gave ori.glnall;'' for the incompressible 
fluid and which was later supplemented by Lin (reference 5 , part III) . 
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At very large Eeynolds nicibers, the relation (2.1) "between the phase 
velocity and the vave length can he considerably simplified. Wlien \ 
is finite and c does not vanish^ \z \ » 1 at large Eeynolds 
numbers. The asyn^totic behavior of the Tietjens function F(z) 
as |zj — is given by (refei'once 5^ part I) 



(3.1) 


and the relation (2.1) becomes 



(3.2) 


where 



is given by equation (2.1h). 


Suppose that a neutral disturbance of nonvanishi.ng wave 

number a_ = ~ and phase velocity c_ > 1 - -1- exists in the 
A s 

° o 

inviscid fluid (lim3.ting case of an lixfinite Reynolds number). The 
phase velocity c is a continuous function of - S, and for a dis~ 
turbance of given wave number tlie value of c at vex*y large 

Eeynolds numbers will differ from c^ by a small increment 2ic. 

Both sides of equation (3-2) can be developed in a Taylor’s series 
in Ac. and an expression for Ac can be obtained as follows: 






-e 



11 + 0(Ac)] 


(3.3) 
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The "boundary condition 



(3A) 


must be satisfied for the inviscid neutral disturbance, and the 
function Ei(ag, Cg, vanishes (eqviation 2.14) . Eecognizing 

that 



reduces equation (3.3) for Ac to the form 


rti/4 

Ac = (3.5) 



From equation (2.l4), 



23 


( 


By equations (2.12) and the boundary condition (3.^), the quantity! 




9Si 


3c 




is evaluated as follows: 


n=l n=l 


2-v(i-sy 


00 


3 






n=0 


r»-l 


ro 


whore the primes now denote differentiation ■'..1th respect to c. Foi' small values of Cg and ttg, the 
quantity • is c^ven .apprordmately by the relation 


wie pi’imeu no^' 

fe) a 




2 - Mq2(i - Cgf 

- ^ef - °sf 




I 

(3.7) 


(3.8) 
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and the expression for Ac is 



Evalv-ation of the integral 



yields the folloving result: 


IV 0/ v'c 


fa /wV 


4., ,^3 


(^c’l 


^7 


[rj 


(In c - irt) + 0(1) (3.10) 


Since the quan 


r<i /v'V 

txty — I r~ 

Vt /_ 


vanishes (reference 8), differ- 


w-c. 


entiation of equation (3.I0) gives 


K '( c . -4 

1 V s- o y , 2 


'o I fd 


w^'Cg^ \8c l^v 




\?V-G 


^In Cg - ljr^4- 0(1) 


vc==c. 


(3.11) 


ThviS; iZ^’^Cg, is approxi-inately real and positive for sttall 

1 

values of c„. With c„ > 1 - — ■, I.P. Ac xauct also "be positive 

o 

(equation (3-9))j therefore, a subsonic distiarbance of vave 
length Xg / 0, idiich is neutral in the inviscid compressible 

fluid, is self-excited in the real compressible fl'iid at very largo 
(but finite) Reynolds numbers. 
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In reference 8 it was proved that a noutral subsonic boundary- 
layer disturbance of nonvanishing phase velocity and finite vavo 
length exists in on inviscid conipressible fluid only if the quan- 
tity ^(p~ ) vanishes for some value of w > 1 - — . If this 

4y \ 3-yy 

condition is satisfied, then self -excited subsonic disturbances 
also exist in the flu3.d, and the Laminar boundary layer is unstable 
in the limiting case of an infinite Reynolds numbor. Sy the exten- 
sion of Heisenberg’s criterion to the coDroresoible fluid, it can be 
seen that, far from stabilizing the flov, the smal3, viscosity in 
the real fluid has, on the contrary, a dest abilizing influence at 
very large Rojnoldn numbers. Thus, any laminar boundary-layor flow 

in a viscovts conductive gas for "viiich the qiacntity ) "vanishes 

dy\dy/ 

for soEie vali,ie of w > 1 - — is \mstablo at sufficiently higli (but 

“o 

finite) Reynolds numbers. 


d f dv\ 

Unless the condition ^ scooe value of 


V>1 


is satisfied, aDl subsonic distisrbances of finite ■wave length are 
damped in the limiting case of infinite Reynolds nxmber, and the 
inviscid flov is stable. Since the effect of viscosity is des- 
tabilizing at very largo Rejn-olds numbers, hovover, a laminar 
boundary flow that is stable in the 3-lmit of infinite Reynolds 
numbor is not necessarily stable at large Reynolds numbers viien the 
viscosity of the fluid is considered. (Soc fig. '^(1) .) In fact, 
for the incomprosoible flvuLd, Lin lias shovn that every laminar 
boundary-layor flow is unstable at sufficiently high Rejnolds 


nuaibers, ■vhothor or not the vorticity gi'adient vanishes (refer - 

dy'2 


once part HI) . In order to settle this question for the com- 
pressible fluid in general terms, the relation (2.1) botwoon the 
complex phase velocity and the wave length at large Ro^Tiolds numbers 

d / dw\ ^ 

must now be studied, for flows in \diicli tlie Quantity* ) does 

dyA^V 


not vanish for any value of w > 1 
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4. aTABIUCCY OF BOUD'IDABY LASER AT LAEGE RESNOIDS NUMBERS 


The neutral subsonic disturbance marks a possible "boundary" 
beti-reen the damped and the self -excited disturbance, that is, 
betveen stable and unstable flow. Thus, the general conditions 
under \dilch self -excited disturbances exi.st in the laminar boundary 
layer at large Reynolds numbers can be determined from a study of 
the behavior of the curve of a against R for the neutral . ^ 

disturbances. IJhen the mean free-stream velocity is subsonic i 


the pliysical sitijation for the subsonic disturbances at large 
Reynolds numbers is quite similar to the analogoiis situation for 
the incompressible fluid. The cui’ve of a a/^inst E for the 
neutral disturbances can be expected to have two distinct asymptotic 
branches that enclose a roj^ion of Instability in the a,E-plane, 
regardless of the local distribution of mean velocity and mean 
temperature across the boundai-y layer. When the mean free-stream 
velocity is supersonic > 1^ the situation is somewhat dif- 
ferent; under certain conditions (soon to bo defined) a neutral 


at any value of the Reynolds number. For this reason, it is more 
convenient to study the case of subsonic and supersonic free-stream 
velocity separately. 


The asymptotic behavior at large Reynolds numbers of the ctirve 
of -a against E for the neutral disturbances is determined by 
the relations (2.19) to (2.22) between a, E, and c for real 
of 0 . For small values of a and c, these relations 
are given approximtoly by 


or a seli'-oxcltod subsonic disturbance 



cannot exist 


a. Subsonic Free -Stream Yclocity (Mq < 1 



u = $y(z) 


(^^. 2 ) 
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(i^.3) 


Vj’c 1 


a = 


Ti u 


- |/l - m/( 1 - c)‘ 


(it.M 


As R — ^co, either z — >■ co or z remains finite tiiile 
"both a and c approach 0. These two possibilities correspond 
to two asymptotic branches of the curve of a against E. 

Lover branch .- If z remains finite as E — >co, then c •— ^0; 
and by equation (4.1), $^(z) — >0. Therefore, z— >2.29 vhile 

u — ->2.29 (fig. 9). From equations (4.3) snd (4.4), along the 
lower branch of the curve of a against R for neutral stability 



(4.5) 


c 


2.29 



(4,5) 


and a — >0 at large Reynolds numbers (fig. 4(1)). 


Upper branch .- Along the upper branch of the curve of a 
against R for neutral stability, z — »oo and 



(4.7) 
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vhile u — )1.0 (fig. 9 and equation (i<-.2)). If the quantity- 

does not vanish for any value of v y 0, then "by equation (4.7) 
c must approach zero as z — . Along 'bhis branchy 






a 


(4.8) 




C te 


a 


(h.9) 


- 




and a — >0 at large Eeynolds numbers (fig. 4(2)). 


v‘ 


On the other hand, if — ( — 

dy\T 


vanishes for some value 


of w = ,Cg >0, then by equation (4.7),. c — >Cg and a — 

as both z and E approach co. Eoif, 



(4.10) 



tions (4.4) and 
agains-t E for 


does not vanish (see appendiac D), then by equa- 

( 4 . 7 ), along the upper branch of -the curve of a 
the neuti’al disturbances, 
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ac' 


(<= - %y 


(^.11) 


Vj^'c 

“ = —-f - - “=)’ 


(i|-.12) 


and c — / 0 ^ a — >a J 0 at large Reynolds nvtafbers . (figs. 4 (k) 
s ^ s 
n>2 “ 


and 4(2)) . If 

by 


dy^ \T /. 


vanishes^ the relation ( 4 . 11 ) is replaced 


1 


E a 



ac' 




(4.13) 


vhich reduces to the relation obtained by Lin in the limiting case 
of an incompressible fluid uhen — ^ 0 , the solid boundary is 

insTolated, and -vt** = 0 for some value of v = c^ > 0 . (See equa- 
tion (12.22) of reference 5 ; part III.) 


d / 

If the quantity — { 

^ V 


3 


vanishes ^ the solid boundary (that 


is, for V =0), it can be shovn from the equations of motion 

. 1 

(appendix D) that — - is al^rays positive - except in the 



2 . 

Lay \T 

limiting case of an incompressible fluid. For small values of y, 


the quantities 


dy \T / 


and ~ are both positive and increasing. 
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For large values of hov/ever, 
w* 

fore — must have a maoimian, or 
T 

of V > 0, and this case is no different from tlie ^general case 
treated in the preceding paragraph. Iii the limiting case of an 

pressihle fluid, vhen V' vanishes at the surface, w ” = — 

since w^’’’ alvays vanishes in this case. Frcm eqtiation (’i-.S) 

the relation between a and E along the upper branch of the 
neu.tral stability curve is therefore 


w' 

T~ 

d 


/ w 


0 , 


dj-V 


physically; there- 

— 1 = 0 for scjue value 
T / 


incom- 





ik.lh) 


which is identical with equation (12.19) in reference part IH. 


Thus, regardless of tiae behavior of the quantity 



regardless of the local distribution of mean velocity and mean 
temperature across the boundary layer - when < 1, the curve 
of a agai:ist R for the neutral distui’bances has two distinct 
branches at largo Reynolds numbers. From plysical coiisiderations, 
all subsonic dist’U’bances must be dairped ivhen the wave length is 
sufficiently sm^ill (a large) or the Reynolds number is sufficiently 
low. Consequently, the two tranches of the curve of a against R 
for the neutral disturbances must Join eventually, and the region 
between them in the ct,R-plans is a region of instability; that is, 
at a given value of the Royno3.ds number, subsonic disturbances Td.th 
wave lengths lying betxreen two critical values and Xg (ct^ 

and are self-excited. Thus, when < 1, any laminar 

boundary-layer flow in a viscous conductive gas is uns tab l e at 
sufficiently hifjbi (but finite) Reynolds numbers. 


The lower brancii of the curve of a against R for the neutral 

disturbances is virtually unaffected by tho distribution of ~- 

dy 

across the boundary layer, but for the upper branch the behavior of 
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the q.miitity 


<iy Vt ) 


Is decisive. 


'V'Jhen 



for some 


value of w = c_ > 0 , the neuti'al suhsonic disturbance passes 
continuously into the characteristic inviscid disturbance c= Cg 
and a = ctg as R — >c». This result is in accordance "VTith the 


results obtained in reference 9 for the inviscid compressible fluid 
and is in agreement "vrith Heisenberg's criterion. In addition, 

all subsonic distxn’bances of finite vave length X > X = — (and 

“ a 
s 

nonvanishing phase velocity 0 < c^ < c^^ are self -excited in the 

limiting case of infinite Reynolds number. On the other hand, 
d / v'\ 

"When ~ I ~ j doss not vanish for any value of vr > 0, then 

except for the ” singular” neutral disturbance of ?,ero phase velocity 
and infinite -wave length (c = 0 and a = O), all disturbances 
are damped in the inviscid compressible fluid. This singular 
neutral disturbance can be regarded as the limiting case of the 
neutral subsonic disturbance in a real coirpressible fluid as R — ^cd 


b. Supersonic Free-Stream Velocity > ij 

V/hen the velocity of the free stream is supersonic, the sub- 
sonic boundary-layer disturbances must satisfy not only the differ- 
ential equations and the boundary conditions of the problem but 

also the physical requirement that c > 1 - ~. The asymptotic 

^^o 

behavior at large Reynolds numbers of the curve of a against R 
for the neutral subsonic disturbances is determined by the approxi- 
mate relations (b.l) to (b.b), vlth the additional restriction 

that c > 1 - As c — >-1 - a — >0 by equation (4.4); 

o Mq 

therefore R — >.a. hy oq^^ation (4.3). The corresponding value (or 
values) of z is determined by eqiratlon ( 4 . 1 ) as folloivs; 


0 ^(z) = v(c) = V 






•lfrC=l 


(4.15) 
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d /w’N 

Now from physical considerations, — ( — ) < 0 for large 

va3-ues of y. Therefore, if — ( = 0 (changes sign) for 

fiy \T / 


dy \T J 

value of w = c„ > 1 - -i— , then, in general. 


_dy Vt yj 

w=l--- 


some 
> 0 


and 


^ \ 


<0 (equation (i4-.15)). From fi;7ure it can he 

■ C=:l 


I't 


N 


seen that in this case there is only one value of z ( zi , Bay j 
corresponding to the value of !^^-(z) given hy equation fh.l5). 

From equations (4.2) to (4,4), along the lower branch of the curve 
of a against B for the neutral distiu’hances, 






\3 


- 

M / 
0/ 


a 


(h.l6) 


a « 


l~ — r 

-. -3 - . 

vi f V 'ty 


(t.17) 


and c 1 — 




at largo Eeynolds numbers (fig. 4(k)). The tippei- 


branch of the ciirve in this case is given by equations (4.11) 


and (4.12), or by equations (4.13) and (4.12) if 


£2_/w' 

\T /± 


vanishes, with c — >c_ >1 and a — 4- 0. 

8 j-r 3 ' 

o 
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If 


-{-) ' 

<3y \T / 


'/■anishes for v = 1 - — , then 2 — >•<» as E — 

M., 


along the \ipper "branch of the cxrrve of a against E for the 


neiitral disturbances, and §^{z)- 


W-, 


Now a • 


0 


B 


1/2“ V 


1 ’ 

as c — >• 1 in this case also (eq,uation (4.17) with u-, = l.O) 

M -L 


so that 


B 


2(v^’) 


pt 2 

o 


o II oh. 


I'h 


^ /!i\l ? a- 


dy‘ 


(9li 


(4.18) 


Along the lower branch of the curve of a against B at large 
Beynolde numbers, a K, and c are connected by equations (4.l6) 
and (4.17), vd-th z-, = 2.29 and = 2.29. In spite of the fact 


that 


dy \ T J 


0 for w = 1 - 




a neutral sonic disturbance 


c = 1 - “ j finite wave length does not e:dLst in the inviscid 


fluid unless K^(c) = 


^CX> 

1 

r T o~ 

Jo 1 

I'h- 

L(^ - c )2 ^ _ 


dy is positive. (See 


section 10 of reference 8.) Calculation sho^ra that K 2 ^(c) is almost 

always negative (equation ( 3 .II)); therefore, in general, the sonic 
dist^Irbance of infinite wave length (a = 0) wi.th constant phase 
across the boundary layer exists only in the inviscid fluid (B — > 00 ). 


If 


f w 




dy \T j 
is certain that 


does not vanish for an;^'’ value of w > 1 - 


M. 


it 




_dy 

V jj 


< 0 and by equation (4.15) 


M, 


o 
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$i(.) >0. 

c=l- — 

1^0 


When V <0.580 (approx,), there are tw 
M 

o 


values of z ^z^ and z^, say, with z^ > correspondins to 

tile value of fj^(z) given by equation (4.15) . (See fig. 9.) Along 


the two asya^jtotic branches of the ciurve of a against E for* the 
neutral disturbances, a, E, and c are connected by relations 
of the form of equations (4.l6) and (4.I7), with z and u replaced 
by Zg end Ug, respectively, along the lower branch and by z^ 

and u^, respectively, along the upper branch . At a given vftlue of 


the Mach number, the value of v is controlled by the thennal condi- 

1-L. 

o 

tions at the solid surface. (See section 6.) "Wlion these conditions are 
such that V s 0.58o_, then z,- = z , and the tw as5Ti5)totic branches 

1-—- "" ^ 

M 


o 

of the curve of a against R for the neutral disturbances coin- 
cide. 1/hon V ^ ^0.580 (approx.), it is impossible for a 

neutral or a self-excited subsonic disturbance to exist in the 
laminar boundary layer of a viscous conductive gas at any value of 
the Reynolds number. In other words, if v ^ ^ 0.580 (approx.). 


the lamlnai’ boundary layer is stable at all values of the Reynolds 
number. (Of coxirse, in any given case, the critical conditions 
beyond which only damped subsonic distitrbances exist can be cal- 
culated more acciu:*ately from the relations (2. 23) and (2,29). 

See section 5 ou minim’mi critical Reynolds number.) 


Tho preceddng conclusion can also be deduced, at least qualita- 
tively, from the resudts of a study of the energy balance for a 
neutral subsonic disturbance in the laminar boundary layer. A 
neutral subsonic distui'bance con exist only when the destabilizing 
effect of viscosity near the solid surface, the dating effect of 
viscosity in the fluid, and the energy transfer betw’’oen mean flow 
and disturbance in tho vicinity of the iimer ’’critical layer” « n 
balance out to give a zero (average) net rate of change of the 
energy of the disturbance. (See Schlichting’s cliscussion for 
inccmpressible fluid in roferenco 4.) In reference 8 it is shoimi 
that the sign and magnitude? of tho phase shift in u* ' through 
the inner ’’critical layer” at w = c is determined by the sign 
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and magnitude of the quantity 


1(f)] 


. Tile corresponding 


w^c 


apparent shear stress T^* = -p* u*’v*', idiic]i is zero for v <; c in the 

inviscid cor^iressible fluid, is given hy the following expression 
for v> c (reference B) , 


/- — \p a 'j^c^ ■ r*i 

T * = p * ( ~ Jt — / — 

C *^0X0 7 p / J \ m / 

(Wc’)^ idy\T /j _ 


^'=C 


(^.19) 


If the quantity 


is negative, the mean flow absorbs 


Xf-C 


[s(f)I, 

fd /w'M 

energy from the disturbance; if \TyJ positive, enersr 

passes from the mean flow to the disturbance. In the real com- 
pressible fluid, the thickness of the irnsr critical layer in idiich 

1 

the viscous forces are impoi'tant is of the order of and 


“0/ 

the phase shift in u*' is actually brou^it about by the effects 

of viscous diffusion (of the Quantity p— -) througli this layer. 

dy/ , , 

As shovn by Prandtl (reference 12), the destabilizing effect 
of viscosity near the- solid suirface is to shift the phase of the 
’’frictional” coiiponent dist’jrbance velocity against 

the phase of the ’’frictionless” or ‘'inviscid” component 


in a thin layer of fluid of thickness of the order of 


]. 


E 

ac — - 

u 


1 


By continuity, the associated normal coiJponent v^.y'*’* is of the 


order of f- — 



- s: u. *’ 
i inv 


a 


V- 


(it -^ras shorn, in part 1 of 
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reference 8, that for large values of ctR the "frictional" 
components of the disturbance also satisfy the, continuity rela- 

. Su’*^ * Sv* ’ ' \ ■ 

txon + =-o in the compressihle fluid.) The corro- 
de* , ' by* 

spending apparent shear stress’ T^* = -Pj_* u**v*' is given by the 
expression 



But from equations (2.11) 


and 




Ti^ “ 



{h.20) 


(4.21) 


(4.22) 


Since the shear stress associated with the destabilizing effect 
of viscosity near the solid surface and the shear stx’ess near the 
critical layer act rougtily throuf^out the same region of the fluid, the 


ratio of the rates of energy transferred 


approximately 


fih- , 

du* . 

T* dy* 

. <iy* 

iJO 


by the two physical processes is 
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T * 

C 


C 


M 


Ej_* 


V 


” 5;’- 



(•1.23) 


1 

2 


'iv(c) I 


vh 6 re 



If the q.uantity 


d_ 


when w = say^ 



is negative and s^rfficiently large 


then the rate at ^hich oner 57 is ahsorhsd "by 


the mean flov near the inner "critical layer” plus the rate at idiich 
the energy of the disturbance is dissipated by viscous action more 
than counterbalances the rate at vhidh energy passes from the mean 
flov to the disturbance because of the destabilizing effect of 
viscosity near the solid surface. Consequentlj'’, a neutral subsonic 
disturbance •vrf.th the phase velocity c ^ c^ does not exist; in 

fact, all subsonic dist’u*bancos for vhich c > c^^ are damped. 

When < 1, there is alvrays a range of values of phase velocity 


0 _< c c^ for whic.h the ratio 


TT* * ^ 

^1 


given by eqaiation {h.22), 


is small enou^ for neutral (and self -excited) subsonic disturbances 
to exist for Reynolds numbers greater than a certain critical value. 
■VHien Hq > 1 , however, because of the physical requirement 


that c > 1 - — 0 . the possibility exists that for certain 


thermal conditions at the solid surface 


the quantity 

L^\TyJ 

' 'Mr.C 


is al'vjays sufficiently large negatively (and tho3?efore 



is 
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sirfficiently large.) eo that only damped suhsonlc dlstijrbances exist 

at all Reynolds numbers. Of course. If ~( —) vanishes for some 

ay\T/ 

value of v^l , it is certain that v(c) <0.580 for some 

- 

range of values of the phase velocity l-~<c<c. In that 

M = = o 

o 

case, neutral and self-excited subsonic disturbances always exist 
for E > ^crjjj^-j^ flow is al^rays unstable at sxifficiently 

hig!a Eeynolds numbers, in accordance id.th Heisenberg’s criterion 
as e:ctended to the compressible flixid (section 2). 

A discussion of the significance of these results is reserved 
for a later section (section 6) in •\daich the boliavlor of the quan- 

^ related directly to the thermal conditions 

at the solid surface and the free-stream Mach number. 


5. CR3HISRI0N FOR THE MUilMUM aRIflCAL EEYHOIDS HmiBER 


The object of the stability analysis is not only to determine 
the general conditions under idiich the laminar boundary layer is 
unstable at sufficiently hi^i Reynolds numbers but also, if possible, 
to obtain some simple criterion for the limit of stability of the 
flow (Mnimum critical Rejuiolds number) in terms of the local 
distribution of mean velocity and mean tomperatia-’e across the 
boundary layer. For plane Couette motion (linear velocity profile) 
and plane Poiseuillo motion (parabolic velocity profile) in an 
incompressible fluid, Synge (reference 13 ) iss able to prove 
rigorously that a minimum critical Eeynolds nimiber actually exists below 
which the flow is stable. His proof applies also to the laminar boundary 
layer in an Incoii^yrossible fluid, ^Tith only a sli^rb modification (refer- 
ence 5, part III) . Such a proof is more difficult to give for the laminar 
oundary layer in a viscous conductive gas^ however, the existence, 
in general, of a minimum critical Reynolds number can bo inferred 
from purely pliysical considerations, k study of the energy balance 
for small disturbances in the laminar boundary layer sho'\i?s that the 
ratio of the rate of viscous dissipation to the rate of energy 
transfer near the critical layer is 1/R for a disturbance of 
given wave length while the energy^ transfer, associated with tlio 
destabilizing action of viscosity near the solid surface bears the 

ratio 1 /)/r to the energy' transfer near the critical layer. Thus, 
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the effects of viscouB c*J.ssipation ■'.dLl?. predomlmto at sufficiently 
low RejTiolds numbers and all subsonic disturbances ^d.ll be damped. 

The two distinct aej-rmptotic branches of the curve of a, against R 
for the neutral disturbances at large Reynolds numbers must Join 
eventually ( section 4) and the flow is stable for all Reynolds 
numbers loss than a certain critical value. 

An estimate of the value of Rf,-v* > t!hich will, serve as a 

^■^min 

stability criterion, is obtained by taMng the phase velocity c 
to have the maximum possible value c^ for a neutral subsonic 

disturbance, that is, for c > c^ all subsonic disturbances are 

damped. This condition is vei’y noai’ly eg.uivalent to the conclition 
that oR be a minimum, which vras employed by Lin fox- the case of 
the inccanpresslble fluid (p. .--85 of reference 5^ pax-'t HI) • The condi- 
tion c = occurs -vdien is a maximum: -fchat is. •'.•Jhen =0.580, 

= 3.22 and = 1.48 (fig. 9). fb-o corresponding value 

of c - Cq can be calculated from the relations (2.19) to (2.22), 

Neglecting terms in (X is usioallj'- very smell) end taking u=l. 50 

gives 

^-^(z) s: j_l - 2X(c)J v(c) (5.1) 


where 


v(c) = - Jt 


wp'c 

't 7* 


d. f w ' 

(w')3 ay\T/ 


(5.2) 


w=c 


;nd 


X(c 




(5.3) 


It is only necessary to plot the quantity'’ (l - 2X)v against c 
I for a .given laminar boundary- layer flow and find the value of e = 

I for which (l - 2X)v = O.58O. The coiTonponding Vvalue of oE is 

determined from the re?.ation 
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(5A) 


and this value of dR is very close to the minirautp. value of ctR. 
A rou^ estimate of the value of a for c = Cq is given by the 

following relation (equation (2.27)): 


~ - ^ 0^(1 - ( 5 . 5 ) 

This estimated value of a is^ in general, too small. The 

following estimate of is obtained by maldns an approximate 

‘^^min 

allowance for this discrepancy and by talcing round numbers: 


“1' 

=0**^ - - %)' 


E 


cr 


min 


( 5 . 6 ) 


or 




E^ 


,76 


cr, 


min 




(5.7) 


For zero pressure gradient, the S-lope of tiie 

/bw\ 

the surface ( J is given very closely by 


velocity profile at 
(appendix B) 
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Wi. 


•B 


•^3. 

0.332 


T 


1 


Therefore 




^^Edn -^1 



[T| 

(“c)j 

|1.76 

( 5 . 8 ) 


1 - 

- = 0 ) 

^2 


The erxprsssion (5.8) ■useful as a rou^ criterion for the dependence 


of E 


on the local dis'tribution of mean Yelocity and mean 


cr. 


min 


temperatvire across the houndarj- layer. It is irnediately evident 

1 


that Eg _ — > CO vhen c^ 


~>1 - 


cr. 


min 


M. 


• . When [(1 - 2X)v] ^0.530, 


c=l-: 




the lonin^r "pgundary layer is stable at al3. values of the Eeynolds 
T ujm'her . fahls condition is an improvement on the stability condi- 
tion V . 0.580 (approx.) stated in section i-.) 

1-L 


In the follo^dng tables and 3 ji figures 5 and 6(a) the estimated 
values of Eg^^ given by equation (5.8) can be compared -with the 


values of Eg 


min 

taken from the calciH.ated curves of Og against Eg 

^^min 


for the neutral disturbances. For the insizlated surface, the values 
are 


U2 
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Wo 

'’^0 

®(°o) 

°^min 
(est.) . 

(fig. 

0 

0.4186 

1.0000 

195 

150 

.50 

. 4^100 

1.0408 

170 

136 

.70 

.4600 

1.0782 

150 

126 

.90 

.4850 

1.1254 

129 

115 

1.10 

.5139 

1.1803 

109 

104 

1.30 

.5^50 

1.2406 

92 

92 


For the noninstaated surface >,iien = O.70, tlio values are 


^1 

0 ' 
0 

1 '(= 0 ) 

°^min 

(est.) 

^^0 

(fig. 4 ) 

0.70 

0.1672 

0.7712 

5377 

5150 

.80 

.2619 

.8716 

1463 

lldA) 

.90 

. 339 ^^ 

.9562 

524 

523 

1.25 

.519^ 

1.1449 

1 

89 

63 


The expression ( 5 . 8 ) for Bq gives the correct order of 

laagnitude and the proper variation of the stability limit with l^ch 
mmiber and with surface tenperature at a given Mach number. 


The form of the criterion for -blie minimum critical Reynolds 
number (equation ( 5 . 3 )) and the results of the detailed stability 
calculations for 8 eve:.’al representative cases (figs. 3 and 4 ) show 
that the distribution of the product of tlie density and the 

vorticity p— across the boundary layer largely determines the 


l.imits of stability of laminar boundary-layer flow. The fact that 
the ^ proper Reynolcls number that appears in the boundary-layer 
staoility calculations is based on the kinematic viscosity at the 
inner critical layer (where the viscous forces are important) 
rather than in the free stresDi also enters the problem,, but it 
amounts only to a numerical and not a qualitative change when the 
usual Reynolds number based on free-stream kinematic viscosity is 
finally computed. Whether tlae value of 2 ^ for a sdven 

°^min 
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laniinar totuidary-layer flow Is larger or smaller tJaan the value 
of Eg for the Blasius flow, for eraE5)le, is determined 

®^min 

entirely hy the distribution of p — across the bomdary layer, 


If the quantity 


d f dw 


;) 


is negative and large near the solid 


tively more stable than the Blasius flow. If the quantity 


surface so that the quantity (l - 2X)v(c) reaches the value O.58O 
when the value of c = is less than 0.4186, the flow is rela- 

iy 

is positive near the solid sxjrface, so that (l - 21)v(c) = O.58O 
wiien w(or c) >0,4l86, the flow is relatively less stable than 
the Blasiiis flow. Thus, the question of the relative influence 
on “ 


(•I) 


E0 


cr, 


of the kinematic viscosity at the inner critical layer 


min 


dv 

and the distribution of p— 

dy 


across the boundary layer, which 


remained open in the concluding discvissions of reference 8, is now 
settled. 


The physical basis for the predominant influence on Eg 

dw 


of the distribution of p 


dy 


‘^^min 

across the boundai-y layer is to be 


found in a study of the energy balance for a subsonic boiuidary-layer 

disturbance (section 4) . The distribution of p — determines the 

djr 

maximum possible value of the phase velocity o^ or the maximum 

possible distance of the inner critical layer from the ^solid surface 
for a neutral subsonic disturbance. The greater the distance of 
the inner critical layer from the solid surface, the greater 
(relatively) the rate of enerQr absorbed by the mean flow from the 
disturbance in the vicinity of the, critical layer (equations (4. 21) 
and (4,22)). l^Then c^ is largo, thei-eforo, the ener®r balance 
for a neutral subsonic disturbance is achieved only when the 
destabilizing action of viscosity near the solid surface is rola- 

1 3/2 

tively large or, in other TOrds, ■vdien x c ' is large 




E. 


Ur 


and the Eeynolds mmiber E , 


which is very nearly equal to E 




is correspondingly small. On the other hand, when c^ is small 
and the inner critical layer is close to the solid surface, the rate 
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at ■v^ich energy is absorbed from the disturbance near the critical 
layer is relatively small and the rate at wMch energy passes to 
the disturbance near the solid surface, which is of the order 

of — .1 — ^ is also relatively small for energy balance; conse- 



quently I^crinin 

6. mSICAL SIGWHTICMCE OF RE3UI.TS OF GTABILm' ANALYSIS 

a. General 


From the results obtained in the present paper and in refer- 
ence 8, it is clear that the stability of the laminai- boundary 
layer in a coii5)ressible fluid is governed by the action of b oth 
viscous and inertia forces. Just as in the case of an incompressible 
fluid, tlie stability problem cannot be understood unless the viscosity 
of the fluid is taken into accomt. Thxis, t'iiether or not a laminar 
boundary- layer flow is ixastable in the inviscid compressible 
fluid (E — ^co), that is, vhether or not the product of the density 

and the vorticity p~ has an extremimi for some value of w> 1 

dy 

there is alxv’ays some value of the Eeynolds number E^,, below 

''^min 

which the effect of viscous dissipation predominates and the flow 
is stable. On the other hand, at very large Eeynolds mmbers the 
influence of viscosity is destabilizing. If tlio ftroe-stream 
velocity is subsonic, any laminar boundary-layer flow is imistable 
at sufficiently hif^ (but finite) Eeynolds numbers, vhether or not 
the flow is stable in the inviscid fltiid ’.iaen oiiLy the inertia 
forces are considered. 


The action of the inertia forces is more decisive for the 
stability of the laminar boundary layer if the fi'oe-stroam velocity 
is supersonic. Because of the physical requirement that the rela- 
tive phase velocity (c - 1) of the bo\mdary-layer distm'bancos 


must be subsonic, it follows that c > 1 > 0 and the quan- 

I-I 

o 


tity 


'W=C 


can be largo onou^ negatively under certain 


conditions so that the stabilizing action of the inertia forces 
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near the inner critical layer (-vdiere v = c > O) is not overcome 
hy the destabilizing action of viscosity near the solid avirface. 

In that case^ mdamped disturbances cannot erist in tlie fluid, and 
the flox'j' is stable at all values of the Seynolds number. 

Eegardless of the froe-stream velocity, the distribixtion of 

d.V“ 

the product of the density and the voi*ticity p~- across the 

fly 

bovindary laj'er determines the actual limit of stability, or the 
minimum critical Reynolds number, for lacn'.nar boundary-layer flow 
in a viscous conductive gas (cq,im,tion ( 5 . 3 )). Since the distri- 
dw 

bution of p— across the boundary layer in turn is determined by 

fly 

the freo-etream Mach munber and the thermal conditionvS at the solid 
surface, the effect of those physical parameters on the stability 
of laminar boundary-layer flow is readily evaluated. 


b. Effect of Free-Stream Mach Wumber and Thermal Conditions at 


Solid Surface on Stability of Laminar Boundary Layer 

The disti’ibution of mean velocity and mean temperatm-e (and 
) across the laminar boxuidaiy layer in a viscous 


therefore of 


conductive gas is strongly influenced by the fact that the viscosity 

of a gas increases with the temperature. (For most gase 3 , |i oc T^ 
(m = 0.76 for air) over a fairly mde temperature range.^ Wien 
heat is transferi’ed to the fluid tlirough the solid suirface, the 
temperature and. viscosity near the surface both decrease along tho 
outward normal, and the f 3 .uid near tho surface is more retai^ded by 
tho viscous shear than the fluid farther out from the surface - as 
compared -vri-th the isothermal Blasius flow. Tho velocity profile 
therefore always possesses a point of inf 3 .ectlon (wiicre w” = O) 


when heat is added to the fluid throuj^i the solid surface, provided 
there is no_ pi'essure gradient in tho direction of the main flow. 

■w*' v*T’ 0.4 4. d / dv\ 

— ~, the auantity — I P } 
t 2 dy \ dyy 

dw 

«iifl- P— - has an extremum at some point in the fluid. On the other 


Since 


d / dw\ _ ^ 

fly \ fl-y/ T 


vanishes 


fly 


hand, if hoat is withdrawi from tho f 3 .uid tlirou^ the solid sur- 

du 


o. cT 

face, — 


end — — are both positive near the sta’face and the 
oy 

fluid near the surface is less retarded than tlie fluid farther 


Sy 
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out - 38 coEpared ■vri.th the Blasitia flow. The velocity profile is 
therefore more conven near the svirface than the Blasius profile. 

As pointed out in section 11 of reference 3, the infD.uence 
of the variable viscosity on the behavior of the pi-oduct of the 

d.v 

density and the vorticity p — can be seen directly from the eq.ua- 

dy 

tions of motion for the mean flow. T-Jhen there is no pressure 
gi-adient in the direction of the main flovr, the fluid acceleration 
vanishes at the solid surface, or 



( 6 . 1 ) 


and 



Thus, when heat is added to the fluid through the solid surface 

is positive, and the velocity profile is concave 


' ")■ (S) . 


near the surface and possesses a. point of inflection for ooxao value 
of w > 0; when heat is withdrawn from the fluid > 0^, ' ^ 


V’'- 

is negative, and the velocity profile is more conver near the siu-faco 
than the Blasius profile. 


The behavior of the quantity ^ 

ex' 


\,T* 5y* 


Vy dy^' \ J 


is 


parallel to that of 


3‘~U'^ ' 


I'-rom equation (6.2), in nondimens ional 


47 


form, 


NACA TN Ko . I 36 O 


[dy \ ay/j^ L^y \T )\^ ^^2 1 1 


( 6 . 3 ) 


Differentiating the dynamic evjuations once and mald.ng use of the 

/w* 

enerv^y eq.uation gives the folloving expression xor 
(appendix D): 


Lay‘-\T/J3_ 


d^ Ar'^ 


ay^\T 



= a(m -'r 1)(7 - 1 )%^ 


-■1 


(^ll- + 2(m + if -rf • (6. >*) 

if ^ if 


Thus, for zero pressure gradient. 


d2 /v‘\| 


_ay^ \T yi 

Now, if the surface is insulated, the quant: 


is alvays positive. 

vanishes, 


hut 






and 


V* 

T 


5oth increase vltli 


distance from the solid sijrface. Since ^0 far from the solid 

svirface, -1- has a maximm and — f 1-') vanishes for some value 
T <3y\T/ 

of w > 0. If heat is added to the fluid throu/Ji the so3.id sur- 
face < 6\, — is already positive at the surface, and 

VI / dy \T / 


since 


dy^ \T 


> 0, the quantity — vanishes at a point 

ay \T / 


in the fluid "rthich is farther from the svirface than for an insulated 
boundary at the same Mach number (figs. 3(a) and (h)). Conse- 
quently, the value of c = c^ for •vdiich the function 
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(1 - 2 \)v(c) = - ltd - 2 X) 


v^'c 



_(v' )3 ay \.T Jj 


reaches the 


value 0.580 is larger than the value for the insulated surface. 

3y equation (5-8). the effect of adding heat to the fluid throu^ 
the solid siuface is to reduce Hq and to destabilize the 

*^^iain 

flow, as compared with the flow over an insulated surface at the 
same Mach number (fig. 6) . 


If heat is •^•d.thdrawn from the fluid through the solid surface, 

is negative. In fact, if the rate of heat 


' > 0 and 


[5(f)], 


transfer is svifficiently large, the quantity 


dy 


(?) 


does not 


vanish within the boundary layer (fig. 3('b)). The value of c = c^^ 

for which the function (l - 2l)v(c) reaches the value O.580 is 
smaller than for an insulated surface at the same Mach number, and 
by equation (5-8), the effect of withdra^Tlng heat from the fluid 
throu^ the solid surface is to increase Sg , and to stabilize 


the flow, as compared with the flow over an ins^tlated surface at 
the same Mach number (fig. 6) . V^ien the velocity of the free stream 
at the **edge” of the botci6.ary layer is supersonic, ttio laminar 
boundary layer is completely stabilized if the rate at •vdiich heat 
is v/ithdrawn through the solid surface roaches or exceeds a critical 
value that depends only on the Mach number, the Ee^violds number, 
and the properties of the gas. The critical rate of heat transfer 


is that for which the quantity — ( — ) is sufficiently la^’ge 

dy \T J 

negatively near the surface (soe equation (6.3)) so that 


(1 - 2>,)v(c) = 0.580 tdicn c = c.^ = 1 - -i— (sections 1)- and 5). 

0 y ^ 

o 

Although detailed stability calculations for supersonic flow over 
a noninsulated svirface have not been carried out, the function 
(1 - £^)v(c) has been con^iuted for noninsulatod sva’facos at 
Mq = 1.30, 1.50, 2.00, 3.00, and 5.00 by a rapid approximate method 

( appendix C ) . The corresponding estimated values of Eo were 

ca...cul.ated from equation (5-8), and in figui'© 7 these values are 
plotted a^inst Tj , the ratio of siurface temperature (deg abs.) 


to free-ctream torperature (deg abs.). At any given Ifech number 
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greater than unity tho -mlue of No increases rspidly 

°^min 

as c — >1 - liien c differs only sli^tly from 1 - , 

o o % 

the stability of the laminar boraadary layer is extremely sensitive 
to thermal conditions at the solid sm-faco. At each value of > 1 

there is a critical value of the temperature ratio ^1^^ ^°r 

tdiich Ep — . If Tt ^ It , tho laminar boundary layer is 
®^min cr 

stable at all Eejaiolds numbers. The difference botvoen tho 
stagnation-temperature ratio and tho critical-sixrface-temperature 
ratio, 'vdiich is related to tho hcat-ti*ansfer coefficient, is plotted 
against I«Iach number in fif^o 3. Under froe-fli^t conditions, for 
Mach numibors groator than some critical Mach number tlaat depends 
larijely on tlie altitude, tho value of Tg - T]_ is "within tho 

order of ma^it-udo of the difference between stagnation tempera ture 
and surface temperature that actually exists because of heat radia- 
tion from the surface {refox’onces lU and 15) . In other words, tho 
critical rate of heat irlthdrawjil from the flu:ul for laminar sta- . 
bility is within the order of magnitude of tho calctilatod rate of 
heat conduction threuf^ the solid surface which balances the heat 
radiated from the surface under equilibrium conditions. Tho calcula- 
tions in appendix E show that this critical Mach number is approxl- ■ 
matelyl at 50^000 feet altitude and approximately 2 at 
100,000 feet altitude. Thus, for M^ > 3 (spprox.) at 50,000 feet 

altitude and > 2 (approx.) at 100,000 feet altitude, the 

laminar boundary-layer flow for thermal o'qui,librium is completely 
stable in tho absenco of on adverse pressure gradient in the free 
stream. 


blien there is actually no heat conduction throuf^ tlie solid 
surface, the limit of stability of the laminar boundary layer 
depends only on the free- stream Mach number, that is, on tlie extent 


of the aorodyn.am3.c heating” of tho order of 





near 


tho solid surface, A good indication of the influence of the froe- 
stroam Mach numboi- on the distribution of p— across tho boundary 

iy 


layer for an insulated siurfaco is obtained frem a rou<^ estimate 


of the location of tho point at whloh 


d / dw\ 

ciy 


reechos a positive 


maximum ( or 


d"^ / dw 

~o P- 
dy^ \ dy , 


vanishes) . Differentiating tlie dynamic 
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equations of mean motion t-^-rice and maWng us© of the energy and 
continuity eq,UEtions yields the folloving result for an insulated 
siirface: 


d8 /w’\' 

b^ 

_dy3 \T /_ 

2 T- 

1 


hr 

mi-2 


(6.5) 


where t r= &ij 


u. 


II’Jq* X* 


From equations (6.4) and (6.5) the value 


of c at wlilch 


d.v2 V J 


. , a. /w‘\ 

vanisheSj or — | — reaches a maximum, 

<3y \T y 


is given rout^ly for air by 


, , Fd^ /VN 


M, 


G ~ 



r *5 
£1/ 








M 

o 


(1 -f 0. 202514^2^ 


1.24 


(6.6) 


. ,, , , b(0.3320) , 

xn vnicli ^2^ » (appendix B), . In other vords^ the point 

**■][ 

in the fluid at which attains a maxiiaita moves farther out 

dy \T / 

from the surface as the Mach number is increased - at least in the 
range 0 (approx.); therefore the value of c for 

1,- > <?- 7*^’\ 

wmcn — ! vanishes and the value of c = c for which 
dv \T / o 

(1 - 2 X)v(c) reaches the value O.580 both increase vdth the Mach 
number (fig. 3 (e)). Sy equation (5.8), the value of Bq for 

°*miri 

the laminar boundary- layer flow over an insulated surface decreases 
as the Mach number increases and the flow is destabilized, as com- 
pared -!d.th the Blasius flow (fig. 5). 


51 


NACA TN Wo. 1360 


c. Results of Detailed StaMlity Calculations for 
Insulated and Noninsulated Surfaces 

From the results of the detailed stahility calculations for 
several representative cases (figs, i)- to 6) ^ a quantitative 
estimate can he made of the effect of freo-stream liach nxmiber 
and thermal conditions at the solid surface on the stability of 
laminar boundary-layer flow. For the insulated surface, the value 


of R 


is Q2 when = I.30 as compared with a value 
<^^min ' 


of 150 for the Blasius flow. For the noninsulated surface 
at M- = 0,70, the value of Eg is 63 vnien Tt = 1.25 (heat 

° ^crmin 

= 126 when Tn = 1.10 (insulated sur- 


added to fluid) , Rq 
face), and Rq 


'cr 


min 


c^min ^ 

= 5150 when Tj = O.70 (heat -v/j-thdrawn from 


fluid). Since R^ =: 2.25 Rq^, (the value of 6 


u * 

o 


V V 


■jdiich 


is proportional to the skin-friction coefficient, differs only 
subtly from the Blasius value of 0.6667) the effect of the thermal 
conditions at the solid surface on R^ is even more pronounced. 

The value of is 60 x 10^ "vAien = O.7O and = 0,70, 

as compared with a value of 5I X 10^ for the Blasius flow 
= 1 and Mq = 0^, For the insulated surface the value 


of R 




cr. 


min 


declines from the Blasius value for = 0 to a 


value of 19 X 10^ at = I.30. The extreme sensitivity of the 

limit of stability of the laminar boundary layer to thermal condi- 
tions at the solid surface "vdion Tj_ < 1 is accounted for by the 

fact that Cq is small when < 1 and Mq < 1 (or is not 

1 


much greater than unity) and Rq « 

'^^min 


(eq.uation (5.8)). 


Small changes in Cq, therefore, produce large changes in Eg 
In addition, -vdien Tq_ < 1, small changes in the thermal conditions 


<^^min 


at the solid surface produce appreciable changes in 
tion (6.3)) and, therefore, in the value of 




(f) 


(equa- 


Not only is the value of Eg^,, affected by the thermal 

^-min 

conditions at the solid sxirfaco and by the froo-straam Mach number 
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but the entire curve of og. against Kg for the neutral dis- 
turbances is also affected. (See figs. 4(k) and 4(7).) Wlien the 
surface is insulated (and / O), or heat is added to the fluid 

(ti = 1.25), Og — >ag ^0 as Pv0 — >co along the iipper branch of 

the curve of neutral stability, in other words, there is a finite 
range of unstable \7ave. D.engths even in the limiting case of an 
infinite Eeynolds number (inviscid fluid). Hovover, a — >0 
as Kq — > 00 for the Blasius flow, or when heat is 'v/ithdram from 

the fluid. This behavior is in complete agreement with the results 
obtained in section 4 and in reference 8, 


A corrparison between the curves of og against Eg for 
Ti - 1.25 and = 0.70 at = 0,70 shows that withdrawing 
heat from tho fluid not only stabilizes the flow by increasing Eg 


but also greatly reduces the range of unstable mvs numbers (ag). 

On tho other hand, the addition of heat to the fluid through the 
solid surface greatly increases the range of unstable *v/ave numbers. 


^^min 


It shouD-d also bo noted that for given va.luos of og, c, • 
and Eg tho t3.me freouoncics of the boundary- layer disturbances 

in the hi^-speed flow of a gas are considerably greater than the 
frequoncies of tho familiar Tollmien x^ves observed in low- speed 
flow. The actual time frequency n* expressed ncndimonsionally 
is as follows: 


n* cocg 

(v)2 

i ou given values of c. Oq , and Kq tho fusg.uency increases as 
the square of the froe-stroam velocity. 


d. Instability of Laminar Boundary Layer and 
Transition to Turbulent Flow 

The value of obtained from the stability analysis 

for a given laminar boiuadary-layor flow is tho value of tho Eeynolds 
mmiber at which self -excited disturbances first appear in the ' 
boundary layer. As Prandtl (reference 12) carefully pointed out. 
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these initial dieturhances are not turbulence, in any sense, but 
slovly groving oscillations. The value of the Eejuiolds numoer 
u-hich boimdary-layor Wstxurbances propagatoO. along bhe surface vd.ll 
be aurplified to a sufficient extent to cause turbulence must be 
larger than Eo in any case; for the insixlated flc^t-plate 

“ crjnin 

flov at low speeds and vd.th no pressure gcadiont, the transition 
Eeyno.lds number found to be tlireo to seven times as 

large as the value of R 0 „v. (references 6 and ,7). The value 

crjjjin 

of 2^0 depends not only on ^crjuin also on the initl.al 

magnitude of the distiu-bances vdLth the most ‘'dangerous*' frequencies 
(those vd.tii greatest amplification) , on the rate of amp].ification 
of these disturbances, and on the physical process (as yet unloiown) 
by vdiich the quasi-stationary laminar flow is finally destroyed 
by tho amplified oscillations. (See, for oxasiple, references 16 
and 17.) The results of tho stability .■’nalysls nevertheless permit 
certain general sti;tements to be made concoming tho effect of 
froe-stroam Mach nu33iber and thermal conditions at the solid surface 
on transition. T’lo basis for those statements is summarized as 
follovrs; 

(1) In many problems of technical Intorcst xn aeronautics the 
level of frec-stream turbilLenco (magnitude of initial disturbances) 

•is sufficiently lov so that the origin of transition is alvrays to 
be found in the instability of the laminar boundary layer. In 
other words, the value of an absolute lower limit for 

C^TD-ln 

transition. 

(2) The effect of the freo-stream Mach number and the thermal 
conditions at the solid surface on tho stability limi.t 

is overvdielming. For exiiSiplo. foi- = O.70, tho value of R0 c r-ny; n 

v?hcn Ti = 0.70 (heat vrf.thdrawn from fluid) is more than 80 times 

as great as the valuo of Eg' when T-i = 1.25 (hoat added to 

c^min 

fluid) . 

(3) The maximv>m rate of amplification of the self-oxcltcd 
boundary -layer disturbances propa^ted along the smface varies 

roughl"^ ss l;\jE0 . (This approximation agrees closely vrith 
‘ / ^ '^^'min 

• tho numcricul results obtained by Pretsch (rofcrenco 18) for tho 
case of an incompressible fluid.) Tho effect of vri.thdrawlng heat 
from the fluid, for cxariple, is not only to increase , and 
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str.'bllize the flow in that maimer hut also to decrease tlie initial 
rate of an^lification of the unstable disturbances. In other i^ords, 
for a given level of fTee-streem turbulence, the interval 
beti/een the first appearance of self"excited disturbances 
and the onset of transition is expected to bo much longer for a 
relatively stable floi^ for which Bq is largo, than for a 

^^min 

relatively unstable flow, I'or \iilch is small and the 

^^min 

initial rate of anplifi cation is large. 

On the basis of these observations, transition is delayed 

increased) by withdrawing heat fnom the fluid tlirough the solid 
surface and is advanced by adding heat to the fluid tlirough the 
solid sun-face, as compared with the Insulated surface at the same 
Mach number. For the insulated surface, transition occurs earlier 
as the Mach number is increased, as compared -vdth the flat-plate 
flow at vei’y low Mach numbers, VJhen the froe-sti*eam velocity at 
the edge of the boundary layer is supersonic, ti-ansition never 
occurs if the rate of heat withdra\ral from the fluid throu^ the 
solid surface reaches or exceeds a critical value that depends 
only on tho Mach number (section 6b and figs. 7 and 8). 


A comparison bet\TOcn tho results of the present analysis and 
moasuroments of transition is possible only when the fcee -stream 
pressure gradient is zero or is held fixed viiilo the froo-stream 
Mach number or the thermal conditions at tho solid surface are 
varied. Liepmaim and File, (reference I9) have measured the move- 
ment of the transition point on a flat plate at a very low freo- 
stream velocity idien heat is applied to the surface. They found 
by means of the hot-i-rire anemomotor that dccl-‘ned 

•^'tr 


from 


5 X 10^ 


for the InGulatod surfoce to a. value of approxi- 


mately 2 x 10^ for 


turbulence 



-■ 1,3^ -idien tho level of free -stream 
was 0.17 percent, or to a value of 3 x 10^ 


when i = 0.05 percent and Ti = 1 , 40 . 

declines from i^•70 (approx.) to 300 (approx.) 
to 365 in tho second. 


Tho value of 
in tho first case and 


'^d McGullou^ (reference 20) observed the variation in 
the transition Eoynoldn number \dien hoot is applied to tho upper 
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surface of an NACA 65 , 2-016 airfoil at the nose section alone, at 
the section Just ahead of the ininiiaum presaiure station, and Tor 
the entire laBlmr run. When heat is applied only to the nose 
section, the ti-ansition Eeynolds number (determined by total-pressure 
tube measurements) -was practically unchanged. Near the nose, 

Ka « the strong f avoratle press^tre gradient in t.x© 

re'^on of the stagnation point stabilizes the laminar boundary layer 
to'such an extent that the addition of heat to the fli^d has only 
a negligible effect. When heat is applied, liovQver, to the sec ..ion 
Just ahead of the minimum pressure point, uiiere the pressure 
gradients are moderate, the transition Seynolds number 


declined to a value of II90 for a?'i » l.l 4 , compared vith a value 

of 1600 for the Insulated surface. When heat is applied to the 
entire laminar run, Ng declined to a value of IO70 ior p » 

XTx 


It -vrould be interesting to invest! ^te experimentally the 
stabilizing effect of a tdthdra^®! of heat from the fluid at super- 
sonic velocities. At any rate, on the basis of tn© results obta.ned 
in the experimental investi^tions of the effect of heating on 
transition at low speeds, the results of the stab?.lity analysis 
give the proper direction of this effect. 


7 . Stability of the Laminar Boundary-Layer Flovr of a Gas mth a 
Pressure Gradient in the Liroction of the Sh'oo Stream 

For the case of an inccanprcssiblo fluid, Pretsch (reference 9 ) 
has shown that even with a pressure gradient in the direction of 
the free stream, the local mean-velocity distribution across the ^ 
boundary 3nyer cou^letoiy dotormines the stability characteristics 
of the local laminar boundary-layer flow at large Reynolds numoers. 
From physical considerations this statement should apply also to 
the compressible fluid, provided oiily the stability of the flow ... 
in the boundary layer is considered and not the possible inter- 
action of the bovjidary layer and the main "external flow. Further 
study is req.uired to settle this (luestion. 

If only the local mean volocity-temporaturo distribution across 
the boundary layer is found to be significant for laminar stability 
in a compressible fluid, the critcrions obtained in the present 
paper and in reference 8 are then iraaediatoly applicable to laminar 
boundary-layer gas flows in which there is a free-stream pressure 
gradient. The quantitative effect of a pressure gradient on laminar 
stability coiad bo readily determined by means of the approximate 
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In a region of snaiJ. or moderate pressttre gi-adients 





say'\ the distribution of p— is sensitive to the thermal conditions 

J djr 

at the solid surface. For example, the chordtJlae position of the 
point of instability of the laminar boundary layer on an airfoil 
u’ith a flat pressure distribution is expected to bo strongly influ- 
enced by heat conduction throu^ the surface . (See reference 20 , ) 

For the insulated surface, the eq.uations of mean motion yield the 
following relation (appendix D) , which does not involve the pressure 
gradient expli c it ly : 


d^' /'dw? 


dy' 


.2 I 


P 

\ WJ 


rt a(m - 1 - 1) (r - l)M, 


1 


2 


° m2 
•^1 


>0 


(T.3) 


The effect of ’'aerodynamic heating” at the siu-face opposes -toe 
effect of a favorable presoui’e gradient so far as tlie distribution 

of across the boundary 3.ayer is concerned (equations (7*2) 

<3y 

and. ( 7 . 3 )). The relative quantitative inf.luance of these two effects 
on laminar stability can only be settled by actual calculations of 
the .laminar boundary-layer flow in a compressible fluid with a free- 
stream }?resaure gradient. A method for the calculation of such 
flows over an insulated surface is given in reference 22. 

Wien the local free-stream velocity at the edge of the boundary 
layer is supersonic, a negative prossxjre gradient can have a decisive 
effect on laminar stabl.lity. The local laminar boundary -.layer flow 
over an instilated surface, for example, is e:tpected to be completely 
stable i»hen the magnitude of the local negative pressure gradient 
reaches or exceeds a critical value that depends only on the local 
14ach number and the properties of the gas. The critical magnitude 

of the pressure gradient is that which makes the quantity — ( p - — | 

cly^ dyy 

sufficiently large negatively near the STarface so that 


. [1 - DX(c)] 


Wt ’c 


T^ d / v*'' 


L(w')3 dy Vi’/_ 


= 0.580 


TfcC 


when c = 1 , 

M 

o 
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(eq.mtlon (5.7) )> in terms of the distribution 

across the boundary layer. Such calcixlations 

(unpublished) have already been carried out by Dr. C. C. Lin of 
Brown University for the incompressible fluid by means of the 
approximate estimate of Ep.* given in i-eference 5> part IH. 

. c^min 

In any events the qualitative effect of a free -stream pressure 

d-vr 

gradient on the local distribution of p — across the boundary 

dy 

layer is evidently the same in a coE5)ressible flxiid as in an incom- 
pressible fluid. If the effect of the local pressure gradient alone 
is considered, the velocity distribution across the bovindary layer 
is ”fiiller” or more convex for accelerated than for uniform flow, 
and conversely, less convex for decelerated flow. Thus, from the 
resu3.ts of the present paper the effect of a nogstive pressure 
gradient on the laminar boundary-layer flow of gas is stabilizing, 
so far as the local moan velocity-temperature distribution is con- 
cerned, while a positive pressure gradient is destabilizing. • For 
the incompressible fluid, this fact is well established by the 
Eaylei^i-Tollmien criterion (reference 3) , the work of Heisenberg 
(reference 1) and Lin (reference 5)> and a mass of detailed cal- 
culations of stability limits from the curves of a against R 
for the neutral disturbances. These calculations were recently 
carried out by several Gonaan investigators for a comprehensive 
series of presswe gradient profiles. (See, for oxampie, rofei'- 
onces 9 and 21.) " ; 


estimate of *E 0 „;„' " 
‘^^mln 

of the qtiantity p— - 
dy 


Some idea of the relative influence on laminar stability of 
the thortnal conditions at the solid surface and the free-stream 
pressuro gradient is obtained from the equations of moan motion. 
At the surface. 





1 


dp* * 

dx* ® dx^ 

(7.1) 

^ 1 



1 



or 


d f dw\ 


duo* 


m +■ 1 ! ' 1 6^= 

rn f f 

m 2 1 i „ m-!-l , „ 

Tt T, u * dx* 

1 o 


(7.2) 
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It has already been shovn in the present paper that \Aien > 3 

(approx.) the laainar boundary-layer flov vith a uniform free-stream 
velocity is con^jletely stable under frse-fli^it conditions vhen the 
solid surface is in thermal eq.iiilibrium, that is, idien the heat 
conducted from the fluid to the surface balances tlie heat radiated 
from the siirface (section 6b). The laminar boundary-layer flow 
for thermal equilibriica should be completaly stable for > Mg, 

say, vhere Mg < 3 if tliere is a netjative pressure gradient in 

the direction of the free stream. Favorable pressure gradients 
exi-St over the fori-rard part of sharp-nosed, airfoils and bodies of 
revolution moving at supers oni.c velocities, and the limits of sta- 
bility /Eg’ \ of the laminax* boundary layer should be cal- 
\ ‘^^min/ 

ciliated in such cases. 


COWCLUSIONS 


From a study of the stability of the laminar boundary layer 
in a compressible fluid, the folloi>7ing conclusions i-rere reached; 

1. In the compressible fluid as in the incompressible fluid, 
the influence of viscosity on the laminar boundax'y- layer floir of 
a gas is destabilizing at very 3.arge Reynolds numbers. If the 
free-stream velocity is subsonic, any laininar boundary- layer flow 
of ^.s is unstable at sufficiently hi^i Reynolds numbers. 


2. Eegerdlass of the freo-stream Mach number, if the product of 

the mean density and the mean vortlcity has an extx’emum ( — 

\ iy 



vanishes } for some value of w > 1 




(where w is the ratio of 


mean velocity component parallel to the surface to the free-stream 
velocity, and where Mq is the free-stream Mach number) the flow 

is unstable at sufficiently hi^ Reynolds numbers. 


3 . The actual limit of stability of lamincir boundarj''-layor flow, 
or the minimum critical Reynolds number E^ , is determined 

c^mi.n 

largely by the distribution of the product of the mean density and 
the mean vorticlty across the boundary layer. An approximate 
estimate of Eg is obtained that servos as a criterion for 

=^min 
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the influence of free-ctream Mach number and thermal conditions at 
the solid surface on laminar stability. For zero pressure gradient, 
this estiioate i^eads as follows: 




'^^ittin 


^ “o'* f - <=o)^ 


where T is the ratio of temperature at a point -vTithin the boinidary 
layer to free -stream temperature, Tq_ is the ratio of temperature 

at the solid surface to the free-stream velocity, and is the 

value of c (the ratio of phase velocity of disturbance to the free- 
stream velocity) for which (l- - 2X)v = O.580. The functions v(c) 
and X(c) are defined as follows: 


v(c) 


nl — 1 c 

if 


T b /x ow'' 
'SwV V 


Vi=G 



C 


where 


nondiiaensional distance from surface 


4, Ozi the basis of the stability criterion in conclusion 3 and 
a study of the eq.uati 0 xzs of mean motion, the effect of adding heat 
to the fluid througli the solid siirface is to reduce Eg and to 

min 
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destabilize the flov, as conrpared ■with the flow over an Insiilated 
surface at the same Mach number. Withdrawing heat throng the 
solid auT'face has exactly the opposite effect. The value of Eg 
^ '^^min 

for the laminar boundary-layer flow over an insulated surface decreases 
as the I-fech number inci’eases, and the flow is destabilized;, as com- 
pared with the Blasius flow at low speeds. 


5. Wlien the free -stream velocity is supersonic, the laminar 
boundary layer is completely stabilized if the rate at which heat 
is v/ithdraT-m from the fluid throu^ the solid surface reaches or 
exceeds a certain critical value. The critical rate of heat transfer. 


for which ™ — yco , is that which makes the quantity ~ { p — 

'^^min • dy dy 

sufficiently large negatively near the surface so that 


[^1 - 2X(c)n v(c) = 0.580 when c = c = 1 - -i-. Calculations for 

^'o 

several supersonic Mach numbers betw^een I.30 and 5.00 show that 
for Mq 5* 3 (approx.; rhe critical rate of heat withdraival for 

laminar stability is within the order of magnitude of the calciilated 
rate of heat conduction tlirou^ the solid surface that balances the 
heat radiated from the surface under froe-fli^it conditions. 

Thus, for Mq >3 (approx.) the laminar boundary-layer flow 

for thermal oquij.ibrium is coaroletoly stable at all Eeynolds numbers 
in the absence of a positive (adverse) pressure gradient in the 
direction of tho free stream. 


6. Detailed calculations of the curves of wave number (inverse 
wave length) against Reynolds number for tlao neutral bo'undary-layor 
disturbances for 10 representative cases of insulated and non- 
insulated surfaces show that also at subsonic speeds the quantitative 
effect on stability of the thermal conditions at tho solid surface 
is very large. For e^smplo, at a Mach number of O.70, tho value 

^®cr , ^1 " fluid), Eg = 126 

when T^ = i.io (insulated surface), and E^ = 5150 -viien T, = 0.70 

/'V. ^ . ‘^^min -*■ 

(boat T-nthdirawTi from fluid). Since R^* s: 2.25Eg2^ the effect 

on R^ Is oven greater. 

‘^^min 


7 . The resxats of the analysis of tho stability of laminar 
boundary-luyer flow? by ohe linearized metJiod of small perturbations 
musu bo applied id.th care to predictions of transition, wdiich is a 
non_inoar phenomenon of a different order. ¥it2i(2raw7ing heat from the 
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fltiid tlirouf^ the Bolid eurface, however, not only increases Eg 


^^2i3.n 


hut decreases the initial rate of amplification of the self -excited 
distiu:hance 0 , which is roug^.ly proportional to 


•j *^*mln 


addition 


of heat to the fluid throu^ the solid surface has the opposite 
effect. Thus, it can he concluded that (a) transition is delayed 

increasedj hy withdrawing heat from the fluid and advanced hy 

adding heat to the fluid through the solid surface, as congjared with 
the insuD.ated surface at the same Mach number, (h) for the insulated 
surface, transition occurs earlier as tho Mach number is increased, 
(c) when the free stream ve3.ocity is supersonic, transition never 
occurs if the rate of heat withdrawal from tho fluid through tho 
solid surface reaches or exceeds the critical value for which 


1^0 


cr. 


CO ■ 


(See conclusion 5-) 


min 


Unlike laminar inatahi3.ity, transition to turbulent flow in 
the boundary layer is not a purely local phenomenon but depends on 
the previous history of the flow. The quantitative effect of thermal 
conditions at the solid surface on transition depends on the existing 
presBixre gradient in tho direction of the free stream, on the part 
of the solid surface to wdiich heat is applied, and so forth, as 
well as on the initial magaitudo of the disturbances (level of free- 
stream tvrb''Jilence) . 

A comparison between conclusion 7(s)^ based on tho results of 
the stability analysis, and experimental investigations of the 
effect of surface heating on transition at lowr speeds ahovs that 
the resu3.ts of tho present paper give the proper direction of this 
effect. 

8 . Tho results of tho present study of laminar stability can 
be extended to inclvide laminar bovaadary-layor flow;s of a gas in 
wiiich tliero is a pressure gradient in the direction of tho free 
stream. Although further study is i-equlrod, it is presumed that 
only the local mean velocity-temperatiu’e distribution determines 
the stability of the local boundary -3.ay or flow;. If that should 
bo the case, the effect of a pressure grodiont on 3.aminar stability 
could bo easily calciiLated tlirou^ its effect on tho local distri- 
bution of the product of mean density and moan vorticity across 
the boundary layer. 

■V/lion the free-stream velocity at the "edge" of tho boundary 
layer is supersonic^ by analogy id.th the stabilizing effect of a 
withdraw:al of heat from the fluid, it is expected that tho laminar 
boundary-layer flow is coDq)lotely stable at all Eeynolds numbers 
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vhen the nerjative (favorable) pressTjre gradient reaches or exceeds 

a certain critical value that depends only on the Mach number and 

the properties of the gas . The lajjiinar boundary-layer flow over a 

surface in thennal e(3.uilibriuiii should be con^iletely stable for 

M > M , say, where Mg <3 if there is a negative pressure 
o s 

gradient in the direction of the free streaia. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., September 5, 19^6 
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APPEIiDIX A 


CALCULATION OF BJTEGRALS APPEAEBB ItT TBE BWISCID SOLUTIONS 


In order to calculate the limits of stability of the laminar 
houndary layer frcsm relations (2.21) to ( 2 . 29 ) hetvoen the values 
of phase velocity, •wave number, and Reynolds number, it is first 
necessary to calciilate the values of the integrals K^, H 2 , 

Kp, and so forth, ■jvhich appear in tho e:g?rossions for 

the inviscid. solutions cp;i(y) and 9p(y) aud their derivatives 

at the edge of tlie hoxmdary layer. These integrals are as follo-vra 
(equations ( 2 . 13 ), (2.9), and (2.10)): 


E^(c) = 



f-r- 

U-^1 






P'2 



V) = ^2^i-^5=i 
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^13(0) x= K^Eg - K- 


fy(v-c)2 

ay I 

(v-c)2 L , T. 


T-M„2 (v-c)2 
(v- c)^ 


and so forth. 

Terrns of hi^or order than a- in the scries eypress"^ ons 
for 9^ end 9^ are no.r^lected, Vhen a < 1, the error involved 


is small becaxise the terns in the series decline lil-cc — B\ 


ven 


for a > 1, howver, this approxlmatioix is Justified, at least for 
the values of c that appear in the stability calcui.aticnc for 
the 10 representative cases selected in the present naner ‘i^*nr 
exairole, the leading term in K.P. 3 (c), ^here- ‘■^='2,3.'". 

if 3 

is approximately _ j ^^Itiplied by the loading term 

in E.P. N^Cc). The quantity in the brackets is at most 0,12 in 
the present calculations; for exair^jle, ' B.p, N^(c) ^ 0,06 E,?. I'I,(c). 

Moreover, E.P. N^^(c) . (l . c) E.P. Sirdlar apprornitV 

relations exist beUTeen E.P. Ii2j,(c) and' E.P. lL(c); and, in 

addition, E.P. M (c) =: (l - a) ~ E.P. ir.(c) r O.Olp E.P. N^(c), 

o o J ^ 

at moat. 

is inaginaa-y part is calciaated 

is^ At tlie end of this appendix, it is shovn that the con- 
tributions of the ilmaginai’y parts of S^, and N^ are 

negligible in comparison vdth the contribution of I.P. ^3,(0). 


General Plan of Calculation 

The method of calciaation adopted must take into accomt the 

fact that the value of ~ / n^' ) -,-i- -4. 

^aiue 01 ~ j at tile point y = y^, Tdiere w c, 
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strongly influences the stability of the laininar homdaiT- layer. 
Accordingly, the integrals are broken into tvo partsj for exan5)le, 


K^^Cc) = 


ny^ 


UYl 


T 


{ 1^2 


(v - c)^ 


dy + 


T 


Jyj 


dy - M, 


- (c) + 


vhere y^ > y^. The integral Kii(c), 

is calculated very accurately, viiereas 
more approximate method as follovs: 


■^diich involves 

Ki2(c) 


“d /dwV 


IfrC 

is calculated by a 


K,,(c) 






(1) 


This int©{5.’al Is evaluated as a power series in c. The 
velocity profile w(y) is approximated by a parabolic arc plujs a 
strai^t-line segaent for purposes of integration. In the more 


complex integrals 

r.y 


E 


Mn 


grals 


Oy 


j 


(w - c) 


2 


2^ 

dy and 


axxd 
012 


N 


w 


3 ^ 

T 


the indefinite inte- 


(w - c) 


— dy are evaluated by 21 


or 4 l point numerical integration by means of Simpson's ruJ.e. The 
values of w(y) are read from the velocity profiles of figures 1 
and. 2 . The value of y^ - y^ = a is 0,40 in the pi’esent series of 

calcifLations; this value is chosen so that the point y = y. is 

never too close to the singularity at y = y^. Take 

ny-1 


KilCc) = 


T 


i/yi (v - c)^- 


;■ ay 


(2) 
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•ly 


The integral -1 (0)^ or the indefinite integral 

,.yi 

tixat appears in Eg^, and is evnlnated by 03q>anding the 

integrand in a Taylor's aeries in y - y^ and then inte.grating the 

aeries 5 emi by terai. The path of inte0:’ation must be tohen below 
the point y = y^ in the comploz y-plane. 


Instead of calculating the values of the velocity and ten^jora- 

turo derivatives and T^^^^ directly, it is siiapler to relate 

these derivatives to their values at the surface by Taylor’s series 
of the form 



The derivatives at the surface and T^^^^ are calcuJLated 

from the equations of moan motion (appendix B) , 

The integral K^^(c). for example,, is fjna3.1y obtained as a 
power serxes in y^ - y^ = a and in y^ - y^ = a - plus terms 

involving log a. The pliase velocity c is related to a by 


c 



3 ‘ 



where 


\ " 


wi 


(k) 




Terms up to the order of 
terms involving 


S' are retained :Ln order to include all 
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Detailed Calculations 

In order to illustrate the method^ the oyaluation of iC^(c) 
is given in sexae detail, as foUovjs: 

(l) Evaluation of ICt(c): 


(a) Define 


Kl(c) 


Jj^l 




Ku(c) 


17 , 


V 



Now 


T 


(v - 



where 




\!/(y) 1 + (j - 


T 

The function — ~ is developed in a Taylor’s sorios aroi\nd tho 
\|;2 

point w = c as follows: 
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vhere 


Then 


%(c) 


and 







2v ’ 
c 




(k + l)v ’ 
c 




^fr-yc) 

6'c')^ ‘Jy^-yc (y-yc)^ 



^ ’W„ 2 



+ 


4 - 



1 

(k)(k + 




(»1 - -^cf] * 

- - (^1 
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vhere 


- J'c K 

(^} - Vj - (>c - 5^1) = = - ” 



<1 = 

“ 


/'t 


The coefficients / 


are expressed in terms of derivatives 

of T and v at 

y = y^ 

as folloira: 

Define 




1 3 /V 

f^(y) = i -~(~-) k>2 

(k - 1) kl 
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(The method adopted for the calciilatlon of from the 

velocity and temperature derivatives and is £iven 

at the end of this appendix.) 

From the expression for K-j^t(c), 

I.P. ^ I.P. K^(c) 

‘ -A'fc) 


= JC 


A(yi) * •' 


J 1 

^ X i 

^ (2^1^ ! 

5- V 


and 


1 . 2 5 


Kin(c) + r— = C_^ C,a -!- C^CT + . . . + 0^.0 

V^'c o 1 c: 5 


I.P. 




ay a - a 


'■o Ai 


ki 


+ . 


•where 


a = V - V, 
- n. *'^1 


. o 


(/l) 


720 




■=k = ^ + — -Th^^ - «k-n A(yi) 


(k + 1)! 


0 < Ic < 5 


(=5 “ °) 


= 


= 4. a f_(y^) + a\(ji) -i- ^\(yj) - 
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(b) Define 


p72 


E^gCc) = 


T 

(v - c)‘ 


<2y 


ni.o 


Jo. 40 


^(j - ’J\) 


CO 


^ £ (k + l)c" 
Z — ic 

k :=0 


where 


-ll-O J 

a.fy - 7 ^ 

Jo . 40 


The velocity profile w(y) ie spprojamated hy a parabolic arc 
in the interval 0.40< y - 4 - 7-, ^ straight lino 

(v = Constant - '*^4.0 interval y^ - 4 ^ 

The value of y, is detortained by iraposing: the condition that the 

5 

area under the parabolic-arc strai.f^it-line sentient .eauals the az’ea 
under the actual velocity profile w(y) in the interval 

O.HO 4 y - y^ 4 parabolic arc w = 2 -!• m^y - y^'^ + n^y 

is determined by the folloT-ring conditions: 

when y = yj^ < 1^ 

w = 1 

w’ = 0 




73 


>1 


NACA TN Wo . 1360 


vhen 7 - 7 ^ aii'i = 0.40, 

■where is road off the velocity profile of figures 1 and 2. 

V <jy 

The valxie of yj^ is chosen so "tlaat -tiie parabolic arc fits the 
velocity curve w(y) closely over the widest possible range. 

For a = 1, 



Therefore 


’'1(1^2-'- J^a)- [(Ti - 1) ~ “o^](ifca * Jk^i) - — ■^k) 


vdiere 


and 


Ic 

Oo.ito 




fll.O 


ty3-yi c'l- 


= - ■ ' (^3 ~ ^x) 


V 


’'^3)1 


Ic 


L 


7h 
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I, is evaluated by approxj.matin^ v(y) by a parabolic arc as 
follows: - • 



~ (k- 1)A 


m -h 2n^y-yj_) 




2lk-l 


1 ^ 3 '^'! 


‘O.IjO 


2k- 3 ^ (-n) 
2k - 2 A 


■-k-1 


■\diere A = - hln . 

As a control in the calculation of the series e 2 q>ression 

OP 

y” a^(k+l)c^ for that, from 

the definition of and 


lim 

}r ^ CO 




and therefore 


1 k 

ifPj S-Tx 

The remainder after K - terms in the series for is given 

approximately by 
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|(N + 1) term] 


c 






T 


(2) Evaluation of H-j^(c): 


The integrand of this integral is free of singularities in the 
region of the complez y-plane hounded hy y = y^ and y = y^; 

therefore H,(c) is evaluated by purely numerical integration. The 
actual procedure employed for the calculation of integrals of this 
type is as follovra: (The integral H]_(c) serves as an illustra- 

tion. ) 

(a) Define 




00 / 


ifhere 
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and 


b = 5\ 


u.-- 




(b) With the approximation that the viscosity varies linearly vlth 
the absolute teaiperature, the velocity v is tho same fimction of 
the nondimensional stream function ^ as in the Blasius flov; that 
is, 


v= v(^)= 


vliere I is defined by tho relation d^ » pv dn (appendix B) . 
From these relations 

d.Ti = ^ H = '^’iB 

since d^ = dri^. Moreover, 



for c = 1. 
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(c) Finally, from the relations given in (h), 


E,(=) . ±1 


\uo 


■>D. 


nh. 


Wg dTJ. 




2c 


+ c 




ii 


u. 


T'diere is the value of 5 ' 


for the Blasiua flov. For 


U 


* T* 


o 


the insiilated surfaces, idiich is some'v^iat arbitrary, f/^as 

chosen as 5-^J vhereas for the noninsulated surfaces, b^ - 6.00. 
(The value of at ng = ^}.60 is 0.99!50j 't-.hen rig = 6.00, 

= 0 . 9975 . The value of b for the insulated surfaces is the 

value of Ti at vhich v = 0 . 9950 ; >iiereas b for the noninsulated 
surfaces is the value of t) for idiich u = 0.99^r5*) The advantage 


nb. 


of this procedirre is that the integrals 






dn^ are calculated 


'.;o 


once and for all- and the value of H-,(c) depends only upon the 
values of b and c. In fact. 


Pb, 


UO 




\ ho* ^‘A V » ”0* 
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Also, 




oO 


\Tq diig - - 1.730 


and 


nb 


b = 


dT] = 


v/0 


Ohr 


Uo 


T d-n. 


'B 


= b^ + 1.73 (Ti - 1) + 0.6667 mJ^ 


= b^ 1.73 


(Ti . 1) - 


7 - 2 

i M ^ 


^7-1 n 
3^7 i. M 


2.3967 ^ 

2 o 


See appendix B. (Incidentally, the last relation shovra the effect" 
of free-stream Mach member and. themial conditions at the solid 
surface on the ’‘thiclcness” of the boundary laj’^er.) 

(3) Evaluation of H^Cc): 


Hg{c) = 




1/71 


(w - c)2 


d.v 


ly 


iyi 


(v - c)^ 


. T 


a.y 


T ^ I ^ (w- c)^ 2 

-ay I — ay.M^- 


f) 7 o p3’ 


Jy^ (v-c)2 Jy^ T 


(v- c}‘ 


dy dy 


^1 


T 
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Define 


H2^(c) 


072 


07 




J7i ~ ifyi 


(■V7 - c)‘ 


T 


- dy 


Hp (c) = 
^2 


72 


ny , .2 

(w - c) 


uy^^ uyi 


T 


dj>- dy 


(a) The integral H22(c) is evalioated by methods similar to those 
already outlined for the evaliiation of H^(c). Thus 


HagCc) = 


(V - 0)2 


Uyi L-'^i 


T 


dy dy 


',72 

'ly 

^y2 

iy 

■i72 

<iy 

pi/’ dy - 2c 

• dy 

pw dy + c 

dy 

yi 1 

yi 

yi u 

'yi i. 

yi ^ 


07 


p fly 


1_ 


p-lB 

T1-!b 

iiJo ■ lb 


'■'b^ 4 Bb - 2= 




T dT>_ 


Uo 


po 

•v?g dTi^ .-j- c I Tt]^ d 
UO 


where 


T = T, 


(Ti . 1 ) . L.^i 

_ f 


The nine integrals in the expression for H22(c) evaluated by 

numerical integration xising Sirpson’s rule. 
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p72 




T 


(J7i 


(w - c)2 




07 




(v - c)2 


- ay 


07 




07 


L/7i 


(v- c)^ 




07 . 


(v.c )2 , T , -(v-cjS, 

dy.j. I d^-- i 1_ dy 


t»7i 


Jy. (v-c)2 Jy^ T 


Define 




i7^ 


T 


07 


»-i 


ay 


(w - c)^ 




T 


ay 


He (c) = 
2 i2 


072 


T P (w - c)2 

dy j i —L~ dy 


Jy 


T 


The integral H2-,p(c) is evaluated, as follows: 


H2^g(c) 


T72 




T 


(v - c)^ 


- a^- 


07 


^1 


(w- c )2 


ay 


T 




Jyj (v-c)2 


dy - 


yi 


T 


ih 


1 72 


T 






072 


U7 


c)‘ 

T 


ay 
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But 


0^2 


Oyj 


(v - c)2 


cly = 


and 


(v - c)2 


I a. ... E,(c) 


.jyi 


so that 


Ho (c) 
^12 


K (c) n^(c) 




(v - c)‘ 


cly 


0^2 


jy 


(w - c)^ 


dy 


Define 


iy2 


P(c) = 


T 


uy-» 


(w - c)^ 


0(y; c) dy 


Ih 


h2 .,) 


).Ah - c) 


G-(ri; c)drj 


IJkCA TN IIo. ,n6o 


vnero 


O'o 


G-(ti; c) = 


(v - c)‘ 


»)t, 


T 


- dn. 


1b 


/fb 




\r 




- 2c 


T 


dn -r c 


i/'n 


li’l 


T 


and 


T = - 


(t. - iV~ M ^ 

V J- 2 ° 


7 - 1 2 2 

V M V 

cl 


The inbegra.! P(c) Is evaluated by nitmerical integration using 
Slnpson’s rule; the required values of are read directly off 
the velocity profiles of fi inures 1 and 2. Fin£;Lly, 


E2j 2(<^) •= %2(«) - P(c) 


The integral Hp^^(o) is evaluated in ejnctly the same vay 
as %q(c) Tdiere 


r ,\2 f K? 


T 


and 


r.r »• 


f f f 
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B-P- [%!<<=)] = (”o=® " + 4o=’* + 

-f- (sCb-^a^ + c.a^ f d..a^ - 2ab - 3a^c - ^ 

\- 1 1 0 0 0^0/ 


2 / 2 


+ a i‘b^a‘“ + c_a~^ - 2abi 


2 


V‘ 


- 3s c. - Ua'^d, - — + 3a c -i 6a d +.10a'"n 
1 lo O o o 


3/2 ^ 2 2 
a'^('b_a' - 2a"Dr> - 3a^c~ + 3ac, •:• 6a d, -i - 'iad - 10a n + 
y3 <^ 21 l'*^o o o3 

2ab2 + 330^ + - ii-ad^ + + 5an^ + 


+ 0 


?2 ^1-- 


+ la 


/a - 0 



. i 

J 

|03 ■ !, 

-T 

— ■■ „ rf 

1 A'C^i) , 

36 

, a - 0 1 





lA’(^i) 


0(^1) 




1 / ^1' 

.. _ [l 



3 ^o(^l^! 


viicre 
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ooj H -P'1 H -P'1 H 


ITACA T?I Ko. 1^60 







A « 
2 


rniTl rri*^mf 

n ^ -^1 -1 

m T T 

•^1 1 1 


A ’ V 

- 2 ! ) 

; 
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V- 


(k) 








V* - \+i' " ■ •''^2'^ 


2 1 k ^6 
2 < k <.5 


•" = A, - Ar/i," - 2A_*A, ' - Aj'.A, 

: k+1 2^ k 2 k 2 k 


2 < k < 4 




« A A A )i r- lA ft! If It t -. I ^ i'V'. 

^ k+1 ■ 2^k •••"'-2 "^k ~^^2^'k '•^2 "^k "^'2 ^^k 


k= 2 


, (a) (m) , (ml) 

^k "*k+l " 


2 ^ k 

0 < m 




7 ■ 


B, 


{ Ksf ^ "*¥3" 


Bn = A 


>8 = 


Bq' - 
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B 


10 


- 


Finally^ 


R.P. H-(c) = E.P. Ho (c) -r Ho (c) - K %[p^(c) 
^ ‘-il ‘-12 

(k) Bvaluation of M^Cc): 


O^r 


M3(c) = 


(v - c)^ 


4y 


i/Ji 




uy 


T 


.(v - c)‘ 


- f 

V ' 


(V - 0 = 


dy 


uyi 


T 


M^Cc) . H/M2g(c) 


\!here 


M3i(c) = 


(ty. 


uyi 


(v - c)2 


pyg 


T 


dy 


0 


T 


(w - c)^ 


dy 


(\y f .2 

(v - c) 


dy 


■V-! 


and 


M^oCc) = 


nyr 


tiyi 


- c)‘ 
T 




^‘^2 py 


dy dy 


uy uyi 
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(a) The integral M22(c) is evalixated in much the same vaj 
as HggCc)) ■tha't is? 


h^3g(c) - 




Wg dTlg 


Oh, 








nt. 


- 2d 


dr-B 


It, 


t% uo 




nt. 


Tv^ diig dvi3 + 




VO i.nj3 OO 


Oh. 


i^B 


TWg dTI.3 dT^^ 


1^0 


+ c 


dilB 


UO iJT)b VO 


■1h. 




Twg dTi2dn^+ 


Hh, 


dii3 


UO 


Oh. 


J’Ib 




OK 


+ 1; 


nd 




uo {|tib uO 


It). 


B 


\ /f'h 


T^ig dT)^ dn^V 

/ 

/ 

/ 


Oh, 


dri3 


,J0 


Tng dTig 


f‘ho 

d% 

UO t. 


"Jh. 


IAB 


r<h. 


T\^ drjg dflg I + c 


% yO 


dT). 


j CO 


B 


•^h 


U^B 


Trig drig 


where h has the same meaning as in the eva 3 .uation of 1122(c) 


and where 


I . T(»g) . T, - j(T^ - 1) - ^;-d] 'ij - ^ 
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The integrals in M^gCc) are evaluated by numerical integration 
using Simpson's rule. Values of Wg are taken from the table in 
appendix B. 

(b) For convenience, the integral M2^(c) is transformed as 
follows: 

M^^(c) = + M^i^Cc) - M3 i^(c) 

•vdiere 


py 


1431(0) = 


(w - c)^ 


tyi 


T 


dy 


F / ^2 

T ^ (w - c)^ ^ 

dy i dy 


ily 


(v - cr J 


yi 


T 




)y. 


^ (V - c )2 


pya 


u^i 


T 


dy 


T 


U 7 J 


-- dy 


py 


o’yi 


(w - c)^ 


dy 


T 




)yr 


y^ 


(w - c)^ 


py 


dy 


T 




(w - c)‘ 


dy 


py / ,2 

(v - c) 


b 


dy 


yi 


T 


It is recognized that 


py. 


i/^i 




ri7r 


T 


V (v - c)^ 


dy 


)y 


yi 


(w - c)"^ 


• di^ = (c) 

^12 


91 


MCA TIT No , 1360 


Therefore 


- %(=) 


By additional transforraations, the follo-wlno eq.uations are obtained: 


M^^^Cc) H^(c) P(c) - Q(c) 


where 


Q(c) = 


iyo 


uyj 


(vr - o)2 


py 


T 


dy 


T 


iiy 


, (w - c )2 


dy 


^^2 , v 2 

(w - c) 


i)y 


<iy 




or 


Q(c) = 


3 b 


00. 4b 


hr - c)" 
T 


An 


dt) 


T 


Hb 


dn 


(w - c)^ 


dn 


(io.4b (v - c)'^ ,jr| 


The integral Q,(c) is evalxiated by nmaerical integration 
using Simpson's rule; P(c) is evaluated in the calculation 
of Ho2(g). 


The integral (c) is obtained in exactly the sane 

K 'n (c) and ^2 ^ ^ (c); that is, •: 


as 
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E.P. (c) = 


^.5 


V I? “ 5 


a 

+ a 2a" 


3 

a Ij. a \ 

— + a — I 

9 3/ 


^ s 2^ 


1 ^ 1 (^ 1 ) 


(a - r) f^(y^) 


In (a - a) 


Finally, 


E.P. M^(c) = E.P. 1133^^(0) + H^(c) 


( 5 ) Evalviation of N^Cc): 


1^3(0) = 


pyp r -I ■ py 

d. r >2 2^ 

m/ 


(w - c)‘ 


I (w - c)^ 
dy } i_ dy 


I'iJc 


I rn 

l.i7l 


07 


T 2 

V 


(w - c)2 


dy 


^3(0) 


f!72 

dy 

Jyi (W - c)2 , 

0^2 

M, 


iy 


y7i 


(v - c)^ ^ 
ay 

T 


07^ 


T 


Jy - c)2 


d.y 


dy 


Jyi - c)2 ,J 


07 


c 


72 - 7) 




a7-i 

•'I 


T 


dy 


p72 ,oy 




|7l J7i 


(w - c)2 


dy 


f)7p 


T 


(w - c)‘ 


dy 


P72 


17 


+ M. 


(v - c)‘ 


(i7i l 7 i 


(yp - y) dy 
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It can "be shown that the second and third integrals are identical; 
therefore 


N 


^(c) = N3 i(c) - 2 M^% 32 (c) + M^‘^N34(c) 




vhere 




^3,(0 


T 


dy 


07 f ,2 
(w - c) 


ihi 


nyo 


071 


T 




?/y 


(v - c)^ 


iy 


><3,(=) 


NjjCc) . 


T 


t/yi 


(w - c)^ 


ay 


, ^ (v - c)2 

(’'S ■ ^0 — 


Fs 

^3,^0) = dy 


(v - c) ^ .. 

— i — (='2 - 

ail 


(a) The integral N3jj^(c) is svaluateci hy nvmerical integration 
3 .n a manner sind-lar to H2^(c), E^{o) , and li')2(c); that is, 




yjr 


07 




7i 


(w - c)^ 


')y. 


T 






07 


on nh 

o I 


^^\(lo do 


Jo “■■ 1 , 



^h 1 

nn 

1h 

3h 

3t| 

1b \ 

2c 

dr. 

p\r dri 

dr, + c^ 

dt] 

pdr 

dri) 

1 

0 ^ 

0 J 


0 t 

0 1 

/ 


9k 


NACA TN No . I36O 



The integrals in N^j^Cc) are evaluated ty numerical integration 

in a manner similar to that used in the evaluation of M22(c), and 

so forth. Most of the integrals vill already have "been evaluated 
in tlie calculation of H22(c), and M^gCc). 

(h) For convenience, K22(c) is "broken dOTzrx as follows; 


(■'172 


T 




n^p(c)= _ 

3^ / \2 




(w - c)^ , 


T 


(Vs - 7 ) <37 


J T 


0^1 


(v- 0)^ . 


07/ x2 
(v- c) 


^1 


T 


(72-y)ay 


+ 


P^2 m f'^/ ^2 

T (w- c) . V 

--a.v) ^2-y}dy 

Jyi 
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Let 


W22(°) " ■’■ 


where 




iy. 


T 


Jy, (V - 0)^ 


<3y 


ny 


yyi 


(w - c)^ . . 

O2 "0^ 


^32^(0) - 


)y. 


T 


i/yj 


r. (w - c)^ 




n 

I «/ o 


(jyi 


m ra 


“ y) dj' 


py 


Noo (c) = 


2 


T 


Jyj 


(w - c)‘ 


ay 


pyr 


ijy 


(v - c)2 

— (y^ - y) dy 


Now, 


pyr 


Nop (c) = 
^ 3 


nn 


V, (w - c)‘ 

^ J 


dy 


, s2 

(w - c) 


i.y 


T 


[(^2 - 5"l) -^X 


Since yr> - y-| ~ 1.0, and 


1)72 


He) 


T 


pyp 


0^2 


(w - c)^ 


dy 


(w - c)‘ 


d;>' 
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It is found that 


R^g^Cc) = P(c) - P^(c) 

•v^aere 


P^(c) = 


pyc 


(jyj 


T 


(v - c)2 


- cly 


^2 t ^2 

(w - c) 


(y - y^,) 


<iy 




Oh 


Jo. hi 


G^(t); c)dT\ 


and 




Gn(riJ c) = I — T] dri - 2c 
I T 
ijn 






w 3 

— T] dtj + c 


n ^T] 


ni 




IT) 


P^(c) is evaluated by numerical integration UB3.ng SiE5)8on's rule. 
Define 


1132 (c) 


yg 


T 


yi 


iy 


'iy 2 




(v - c) 


- [(yg - yi) - (y - yi)]«y 


Since 


M2 


T 


I I'r^ 

Ivj 


(w - c)2 


dy = Kj_2(c) 
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and 


(» - C)= 


dy = 3i(c) 




it is recognized that 




ny, 


H^(c) - 






^yp/ \2 
“^(v?- c) 


(y- y3^)<iy= 


1 


f)^ 


p(w - c)%Tl 


Oo 


03 T 



1 


2^ f 

p(w^-2cw+c /dri 

(jc b 

00 UO 


0^2 


(w- c) 




^Fl 


^(y-yi)ay= 


P , 

■b‘"\Jo 


o 


T drijj - 2c 
DO UO 


p 3 

T dn^ 


•H C 


1 % 

Uo 




iij. 


■X dT). 


'B 


UO 


The integral 


f^y 2 2 

^ — ' — ^y - y^^ is evaluated by nurierical 

uyi 


integration in exactly the same vay as W^, (c). 

y'f 
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The integral (c) is transformed as folloiv’s: 


= j 2 

~ ■ Jy. 


>yj A7 

" e, (V . O' 


c'yi 


T 


[(^2 - =^0) - - 5 'o)] 


But 




T 


Jyi - c)^ 


dy 


07 


U 7 l 


(» - 0)2 


and 



so that 


1132^(0) . (1 - cf)H2_^^(c) - J2ii^'=) 


vhere 


P 7 


% (c) 




Sy r !ZJ_£I® 


t-'l 


_J_ ay I iU_W. A. . 5 

(v.-c)2 Jy^ T ^ 


i'o) ^ 


The integral evaluated in the same tray as Kj^j^(c). Thus 
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E.P. 


Jp (c) 


1 


a - CT 






a 


5 



^11 

(yi)^o'i 

[yi) 

j 

Bo(^i); 

r 





+ e'*(5S3o + 


/fl., 

V'* 


r 6a^C - lOa^D^ - iia^C. ) + 0 ^ 


- - haC, 


+ lOa^ 



+ c 


k 



a ^ 2 (^ 1 ) 



•v^iere 


!? 

^l(yl)^b*(^l) 

[foC^oT J 

<lo if- ^1(^1) 1 ^2(^1) 

p ^ H- p + 

° « 25f„(y.) ° ^toQ'i) 


. i_ r^’i'C^i) . 

5 " 15 K(yi) 
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Finally 


E.P. 1132(0) = R.F. 1132^(0) + ^32^(0) - N32^(c) 


(c) Defina 






T 


i -1 - =1® 


ay 


py , .0 

(v - c)^- 
i i_ dy 

T 

Jyi 


3 y 2 


T 


Jy 


(v - c)^ 




After several transformations, the ints-gra.! No-i(c) is brougiit 
into G more convenient form 

^31(0) - 1^31^ (°) I ^^11^*^^ ^ ~ ^^313^^^ 


vhere 




N3: (c) = 




T 


iy 5 'l 


(w - c )2 


'• dy 




(;yi 


(v - c)2 


dy 


■■'y. 


iJy 


(w - c)^ 


dy. 


Hot (c) 

al 3 


oyg 




T 


(v - c)' 




py 


i)y- 


(v - c)‘ 


py 




T 


- dy 


y. (v - c)‘ 


The Integral N^n (c) is evaluated, by numerical integration 
^3 

using Sim.paon's rule. Some of the integrations have alread.y beon 
performed. Tire integral is given as 


-N 3 i,(o) =- 


■b' (.jo 


ni> 


'p 


.Vo (" - 


dr) 


’.11 


Or - c)2 


dr] 


Jo. 4 b 


Vl 


yo 


, 4 b Or - c)‘ 


drj 
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The integral is evaluated in ejccctly the same \ray 

as K-n (c) ; that is, 


1 


[£ 


si , ^5 


■=) = |- ^ foC^l) ■ I'ofc(yz) -f-o’ (''l) j * [(h ■ 




2 ln(a - g) 
(a - 0 ) 


A 


/ 


a 


*,A In 0 + 

a - o 6 


(- 1 ) - e. 3'1 P/,(yi) + I n(7i) 


0 < 


f » 

6 o 


Ai) - T -oC^i) * VoC^i) + ^ 


k 


— E - a% 


" “ ‘■oC^i) " t; C(y.t) - =>11 * VoO'i) - 7 


l 3 


, , 3|.^ J,j, „3 . }_ i f^(y^)- 1- f;'(y^)- a(E + y)_ 

= “'[(Pl - «o)fo(»'l) ^o^o'i^l) - ?^o"(=^l)]) 


+ In (a - 0 ) 


B In 0 -f 


s! f f, v£jp f /-y N.i U 

6 /l(-l) 3 fo^uA/ 3 (* 1, 


+ 0 — f 1 




102 


WACA TN No . 1360 


" 'yf A(^i) " A” (^x) - 5 A ' W " - » iviCV; ^ J 


3 2 \ 3 

+ — a C ) + a"^ 
2 


v.W-ivW 


3 ^oAl) ® 


-rA”(A)-^'“*"' 


- cr 


^xi^i) - PoA'(^i) ■" (*o - A>'i(»'i) 


- i A(?i)-f IfoC^'O+'-A'Cyi)]*!- I\(yi) in c 

■ ‘‘AWi’o] ■'■ 1^- A(yi) + § a(j'i>c 


„ 

12 

_2 


1 


^ r \ i (a - a)^cr 

A(^i)io j " —7— 


af 

If 


L3 

fl(A)]^ 

A( 

.^ 1 ) 


(A) 


1 . .A, 

3 


i'(A) ' '•I’l'oAl) - Vo’(^l) I 


- — (a - <j)f 2 (Vi) - -- (a - ct)^ In (a - 


(a - <j)^at^ '(jl) (a - ct) + (a - a)^ In (a - a) 


|i &i(^i)]‘ 
P A(li) 


- 3 A(a>o 


1 / ^ 
1 / (J 


a - a \ 2 o 


(A)-- 


A’(i'i) h 1 / N 

~i TAAi>o AAi) 
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P 

"■o 1 




(" 1 '^ 











3Pl(yijPl-^Po'’l 







r^l 


u 

1 

kr 

(^■1) 

}\ 


where 


A = 



If 


+ 0 


1 

? 


■(^ 1 ) - 




Vi 




\ o 
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®oA’(^i) * %A(5'i) " V A'C=^i) 


Vi(a) + Vo(>'i) - jA(a) 

Evaluation of fij 

The fvinctionB which appear repeatedly in the 

evaluation of the integrals K]_(c), H2(c), and so forth, are 

evaluated in teniis of and T^ (ic) 

as follows: 
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,V'')2(k. l)k: 

\ ^ J 



k1 


fv'fyi') 


a2^ 


^o‘“V.vi) - -(%),, = - k‘“’ «o * " 


(n-2) 

1 


n 


(m - r) Irl 


m (m-r) , ^ 

J. T , T • ' ' ’ - 1 t-l 


1% 


0 < m < 6 


vhere 


% - 

("l')' 

g. ^ 

-2g Ap 

i 


% 

-2 /g .Ap 
\ 0 *2. 


^1 2j 




ml , (m-r-1) . 

g__Ap + . . . •^.- 


(m - r).'r.' ^ ‘^' 
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0 < m < 3 
•^1 

0^ m < 2 
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r, /_L rn ^ J f » 'S + Jl t ’’S, 


36 ooC’^i’'''®°"' 3Vso}+ eo(-i;;|; ' 180 ‘i '3 ' 2ho '1 


1 1 
+ T, 'S„ - T 

600 1 5 3600 


xh) 


•where 


(k) (k) (k) 1 , (k) 

Sj . 3 B 3 - SCg - - 


0 < k < 3 


s, :r _ 1- 3 B, - 12c 4- — -V 2B 0 < ]a^2. 

4 ^5 4 32^ 7 =i_ 


S 4 4- >530^^^ - 15 c, 

P 3 D 5*0 4 


- 5OB 

6 '^3 2 2 


0 !lk ^1 
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% = - 


+ % * ^10 - + r ^ 


Itc 315 

— C^ + 15c -D^ - 225 E, + — Fo 
31 5 "^3 i. 2 


and are defined ae previously. 


(k) (V) 1 ^ (k+ 1 ) 

Og =: - 2 ^-^2 


5 - ^ \^4 ^ ^ ' ^3 / 


1 4 . -D (k) -n (1^+1) 


V = o V ' -3 


4 


(k) 1 4 . (k) ,. 3 „(k) _ 3 (k-K.LV 


5 


2 




0 (■') = B (-•') - 
6 a 2 


^7 = 4 


^8 " V3^-4 


(k) (k) 1 (k+1) 

> =* ^3 ■ J ^2 


+ 75 D 4 
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::: + C, ' 

^ 3 V ^ 


°U ■ '^s J 


D 


,<’^> c<'=> 

!v 3 y " ^ ^ / 


(k) _ 1 /C, (k) (k) _ ^ (i:+l)N 

^5 - -7 6 J 


- D . Id 

"■•2 " ^3 1^ 2 


E 




'* \ 5 


(k) _ (k) 1 (IK-I) 

- ®3 "5 2 


Order of Magnitude of Inoginary Parts of Interrala Hp^ M- ; and 

J 0 

In the detailed atahility calculations the contrioutions of 
the imaginary parts of the integrals and so forth; 

to the function v(c) are considered to be neglligLtle in cca^arison 
■with -blie contribution of "the inaginf<ry ;^rt of K^(c), A calcui’-atlon 

of the orders of ma ‘jiitude of I.P, H 2 (c), I.P. and I.P. 

from the general expressions given in the preceding pages sliovns that 
this step is Justified, at least for ■tiie values of phase velocity c 
that appear in the stability calcifLations. 

For example. 


I.P. HgCc) = I.P. H2 ^^(c) «^‘(^c')\(yc) 


]1C 


NACA TH Nc. I360 


vhere 


'TheTefor© 


A . . A. _A_ * 0(0") . . . 


I.P. H^(c) » . ■)_! 

A. 3 -V «/ 




Til© contri'bixtioii of I.P. IIo(c) to v(c) is apprcxiinatsly equal 

^ 3 ■’ 

rf. r~^ 

to V 


a c' 

3 T f'-r 

_ Iv 1 y. 


Vn 'C 


, vhere — x.p. Kj^(c). The quantity in 


the brackets is of the order of O.O3, at most, in the caloxaations 
of the present paper. (In the opprox'Lmate calctilations of P»p 

®^Bin 

for Mach nunbers very aucii greater than iraity, c becomes large 
1 


because c > i 


M, 


hovover. a, is small vhon c is not much 


greater than 1 and the results of 'the calculations of Pg 


M 


cr, 


min 




based on tlie approximation v(c) I.P. K^(c) are qualitatively 

^1 


correct (fig. 7).) 


From the expression for Ho(c), I.P. II,(c) 




I.P. K,i(0); 


so that the contribution of I.P. N^(c) to v(c) is approximately 


equa], to v 


of 0.06 at tiae most. 


9 2 
a c 


. The quantity in brackets is of the order 
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The iniegixiary part of M^Cc) is considerably omallor. In 

fact, 


6 

I.P. M_(c) = I.P. K-,(c) 


and the contribution of I.P. M,(c) to v(c) is approrLujately 


equal to 



The quantity in brackets is of the 


or6.er of 0.001 at mazimum c. 
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APPETJDIX B 


CAJXraATION OF MSAN-VELOCin Al® MSAII-TEJ.IPERATI1RE DIgn?IBOTION 


ACE0S3 BOIF3BAEY lAYER AJD TBF VELOCIKr AJID TBffiERATURE 


DERIVATIVES AT THE SOI.ID SURFACE 


The nsan-voloclty and Esean-temperattjre profiles for the several 
representative cases of Insulated and noninsulated surfaces are 
calcu-lated bj 3 rapid approxiinate method that givos the slope of 
the velocity profiles et the sm-face vdth a majdLHimi eiror of about 
4 percent in the extreme case, for ■which T^ - 0,J0 and Mq - 0 , 70 . 

The surface values of the higher velocity derivatives and the 
te 25 >erature derivatives req.u3.red in the stability calculations are 
obtained directly from 'the equations of mean motion in terms of the 
calculated value of t.he slope of the velocity profile. The Prandtl 
nimiber is taken as unity. 


Mean Velocity-Tenrperat-are Distribution across Boundary Layer 

In a seminar held at the California Institute of Technology 
in 1942 , the present author has sho'wn that a good first approxi- 
mation to the moan velocity distribution across tho boinadai-y .layer 
is obtained by assuming tliat tlie viscos3.ty varies linearly i.dth 
the absolute teE 5 >erat luro . With this assimiption, the velocity w(^) is 

the same function of the nondimencional stream function ^ . -X 

as in the Blasius case, and the correopondin.g distance from tho 


Actually, the approximation = *rg(^) is the first stage of an 

iteration process applied to the differential equations of mean 


parameter equal to 0.24 for air), and v(t} - 'WqH) + + • . • 

Calculation of for T^ = I .50 and T^ = 2.00 for M — >0 

shoved that tho iteration process is rapidly con*'’ergent; the con- 
tribution of tho second term to tho slope of tlie velocity profile 



is obtained by a s:Irrple quadrature •vdion a=l. 
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at tho sm-face is 5 percent for T^=1.50 and 8 percent for Tt_=2.00, 
In the present calculations the maximum error in the slope introauced 
hy talcing v(^) = is about percent in the extreme case. 

(See reference 15, in vhich the authors make use of a linear 
viscosity 'temperature relation. See also reference 2^.} 


That w(0 = '^(t) for a linear variation of viscosity "'.ri.th 

absolute temperature is seen directly from the equations of mean 
motion in the laminar boundary layor. The equation of cont_.nuity 
is 8.uto]aaticall;/‘ satisfied ty taking 


and 



p* _ 


The stream function 4^* and the distance along the surface t - are 
selected as independent variables folloving the procedure of van Msos 
p-nri the dynamic equation of moan motion becomes for zero pressure 
gradient 




,2 du* 






— 


A 

“• ““ Wj 


Define the nondimensional stream fiuiction ^ by tiac relation 
\J/* 

f , The dynamic equation takes the folloving form; 


Up * u * 3* 

)| o o 

^ dv d / dv'N 

Since 0 1 in the boundary layer, if [i -- T, the dyniamic equation 

T 

this form is idontical "VJith tha ©g.uation for the isoth63rmal B la si us 


J 
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flow, that is, w(t) s or the valuie of the velocity ratio w 

ia eqtial to the Blasius value at the same va3.ue of The corre- 

spending value of t) = y* .1 ^ -fche nondimensional ** distance” 

from the surface, is obtained as follows; 


or 





u * u * r-'- 


'A 

07 


pw ~ 


A 

dT) 



(.0 


A 

pv 



liO 


V 


If cr = 1, the energy end dj-namic equations have a unique inte;p:'al and 




2l 

J 


\r 



M 

o 


os. shown by Crocco. Therefore, 



115 


MCA TN Ko. 1360 


The integrals 


f 


!o 




Me 


Mb 


Vj,'" dri-g are , given in the 


i'C 


following table, and the mean-velocity and mean-teurperature pro- 
files can be calculated rapidly by this method. (The values 

Bw’ 


of 




are used in the approximate calculation of Eg 


cr, 


B 


min 


(appendix C) . ) 




t = 

1 

[ 

"'E 

0 

1 ^^B 

Jo 

©. 

0.00 

0.0000 

0.0000 

0.0000 

0.3320 

.20 

.0664 


.0066 

.0003 

.3319 

.40 

.1328 


.0265 

.0024 

.3314 

.60 

.1989 


.0599 

.0081 

.3300 

.80 

.2647 


.1065 

.0189 

.3274 

1.00 

.3298 


.1660 

.0367 

.3230 

1.20 

.3938 


.2385 

.0630 

.3165 

1.40 

.4563 


.3236 

.0995 

.3079 

1.60 

.5168 


.4210 

.1468 

.2967 

1.80 

.5748 


.5302 

.2064 

.2825 

2.00 

.6298 


.6508 

.2792 

.2663 

2.20 

.6813 


.7821 

.3854 

.2483 

2.40 

.7290 

1 

.9231 

.4648 

.2280 

2,60 

.7725 

1 

.0733 

.5776 

.2064 

2.80 

.8115 

1.2319 

.7034 

.1835 

3.00 

.8460 


.3973 

.8411 

.1618 

3.20 

.8761 

1.5702 

.969'r 

.l 4 o 8 

3.40 

.9018 

1.7480 

1,1478 

.1180 

3.60 

.9233 

1.9306 

1.3145 

.0986 

3.80 

.9411 

2 

.1171 

1.4684 

.0805 

4.00 

.9555 

2.3067 

1,6682 

.0640 

4.40 

.9759 

2 

.6933 

2.0419 


4.80 

.9873 

3.0863 

2.4280 


5.20 

.9942 

3 

.4828 

2.8211 


5.60 

.9975 

3.8812 

3.2180 


6.00 

.9990 

4.2805 

\ 

1 

3.6167 
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With the approximation that n varies linearly mth the 
ahsolute temperature; the elope of' 'bhe velocity profile at the 
solid surface is simply related to -the slope o'f the Blasius pro- 
file. Thiis 


Since \:(t) = 


and 


/3»A 0.332 

W, ' 17“ 


Sv dv dt dw 

ir) “ dt ^ d| 


Sw ^ c)w\ 



dw 0,332 

dy ~ 1 ~ 


where b is the value of rj at the "edge” of "the boundary layer 
(viien V reaches an arbitrarily prescr.ibed value close to unity). 

It is seen that the shear stress at the surface (or the skin friction) 
has the same value as in the .Blasius case 









^1 V 


5n 

Sy 
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The reliability of this approsimation can be judged from the calcula- 
tions of the skin-friction coefficient in reference 2 k, in ^diich 

1^ cc T*^’ From figiure 2 of reforonce 2 k, the value of the skin- 
friction coefficient for an insulated svirface at a Mach number 
of 3,0 = 2.823] is only 12 pei-oont lower than the Blasius 

value and only 2 percent low'er at a liach mmiber of 2.0 {^'^2. ~ 

For the noninsulated surface^ with = 0.25, the value of the 

skin-fi'*iction coefficient at = 0 is only 7 percent greater than 

tlie Blasius valite and 12 percent greater at a Mach number of 3*00. 

Since the shear stress at the sinrface is uncjianged in first 
approximation, 'yio bound’rj'-layer momentum 'lirickness has the same 
value as for the Blasius flow 


e 



0.6667 


The expression for the displacement thickness b-:;- gives a measure 
of the effect of the theriir:! conditions at the solid surface ajtid 
the free-strerm Mach numbex’ on the thickness of the boundary layei’. 
By definition, 


V 


(ico 




u * xf 

O 4 / 


(1 - pw) dri 


0 


S*, 


From the relation between dr) and drig 

Oc: 

[(T - 1) 4. (1 . d,^ 




o 


UO 


= 1.73 V (t^ - 1)1.73 + lA^io.eeej) 


= 1.73 T, 


7-1 


(0.6667) 
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For the Blasiiis flow 



The ’’thickness” of the boundary layer b is given by 


b 5.60 f (t 3_ - ljl .73 + m/(0.666t) 


S“- 

and the form parameter H = — is 

0 


H =-• 2.50 T, + 2 _l_i M 2 
X 2 o 


For the insulated surface. 


= 2.50 + 3.50 M ^ 

\ 2 ° /■ 


Calci'l.ation of Moan-Ve 3 .ocity and Mean -Temperature Derivatives 

Because of the sensitivity of tlie stability characteristics of 

the laminar boundary layer to the behavior of the Quantity ~ fp — 

\ < 3 y/ 

the values of the req.uired velocity and temperature derivatives at 
the surface are calculated directly from the equations of mean 

motion, t/ith p = T (m = O.76 for air). Nc-vf — = f = bpw so 

t ^ 

that the dinaanlc eqmtion is -b — v' = v’j . Since t(0) = t’(0) = 0 
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V '• ^ 1 'V. ' 

J. m , . X J. 

■^1 *■ 


•Where T, ' = v ' {'LlJ:. M ^ - fT, - 1^1, if CJ = 1. In other 
words, the val-ae of v^" is readily computed from the value 
of w^’. In general, is determined from -tlie relation 

- X (Ev')’"'® = (1“''')’'-" 


or 


V. 


(k) 


. X) -i: V. ^ 


. n” 


+ ' -f 

(p 1 -m X 


(k-3).'2: 


# « • 


(h- 1) 


. (i"X 


(s) 


(k-s) 

' -v ... + 

(k - 1 - b)sI ^ m 1 


Hi 


A ( k-1) 


T 


m 


2P “1-- 




(k-2): ,(k-2-r) (r+1) (k-3)(k-2) ^ „ 

r V + . . . + V. 


(k - 2 - r) ]rl 


1 a 


(k-3) 

_ 
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vhere 

•-- b I '.r !- (p - l)p3^ w 


(p - 1) i (p-a-1) . (q) 

:P-r 


(p - a - 1)1<11 1 


+ 


+ 


p = 1 , 2 , . . .5 


and 


P-. == 




Ti- 


T.. 


P " = 




T 3 tJ 

1 J- 


ar. 




m r I f 


^ m2 

-1 


(t”)'' 


_ r: m(m - 1) 


(Sl')2 Ti" 


+ m 


n 


T, 


T. 


(p-1) 


12.1 
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(^y ' ' 


T “ 

n* 


(m, 2 ) + 3 m(m - 1 ) 


Tl'Ti" 


+ m 


T 3 _'" 


Tn 


T-, 


(t”) 


iv 






(Ti')%i'’ 1 ) 


T 


m 


= (n, 3 ) — — + 6 (m, 2 ) — + 


T, 


ni 3 If ^ 

■^1 


[. 


4 T^'T^"’ 


-1 T 

3 (T{'f U m ^ 
^ ' J T, 


If 5 


m 


T,^ T 3 


15(Tj^‘)(Ti")2 


-h 10 


')\ ' ■ •] + [lOT^"T^” ' + 5I1 'tA’'] + .» ^ 


(ipm) - 

^Z„. = {m, 5 ) 

T. 


D 


15 (:n, M [■‘ 5 (T,')"(V)‘ 


m ^ m 

■^1 -^1 


20(t^ ' ' '] f I 60 T^ 'T^”T^ ’ ' ’ + 15(T3_ ’) + 15(Ti")^] 

.J rn ^ ^ "■ -< 


m(m - 1) 


+ __ — _ — : [ 10/T ' » ’')2 + j _ 5 T '*1 i 

^2 I \1 J ^il 


»if •vl 

“1 J 


+ & 3 _'Tn ■ I + m - 


■^1 

T 


Vi 


1 
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- 5(7 - 


VI 


= av,^'^ - (7 - 1)M^ 


[l0(vi>”)2 


+ 15w 


"t7 

1 


6w^ ' 


Wn 


'] 


Each ve-locity derivative is determined frco the knovledge of all the 
preceding derivatives. 


I 


124 


riACA TN No . 1360 


(m, 2) = m(m - l)(m ~ 2) 

(m, 11) = m(m - l)(n - 2 ) . . . (m - n) 

El = 0.76 
(m, 1) = -0.1824 
(m, 2) = 0.226176 
(m, 3) = -0.506634 
(m, 4) = 1.641455 

( m , 5 ) = -6.959939 ■ 

where 




T 


= aw;j_'” - 3(7 - DM^^w^'w^^” 


^3. ®'^1 


- (7 - i ) k / [3(^1")^ + 


Wn 
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APPENDIX C 

PAPID APPROXIMATION TO THE FHIOTION (1 - 2 X)v(c) AND THE 

CRITICAL REYNOLDS ITUMEER 


In section 5 ; a criterion vas derived for tlie dependence of 
tlie lainimum critical Reynolds nuntei- Rg^^ on the local distri- 

bution of mean velocity and mean temperature across the boundary 
layer. It %-ras found that 



vhere c^ is the value of c for \ 4 iich (l - 2 X)v(c) = O.580 and 



1 


c 
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A rapid method for the calculation of the function (l - 2X,)v(c) 
and the minimum critical Reynolds number is developed by maMng 
uee of the approximation that the viscosity varies linearly vrith 
the absolute teurperaturo (appendix B) . (Since tho effect of 
variable viscosity on the mean- velocity profile is overestimated 
in this approximation, the values of Eg (fis* 6 (a)) calcu- 

°^m;ln 

latod by this method are lower than the values calculated for p = 
when heat is added to the fluid throuj^ the solid surface and higher 
■!?hen heat is '^dthdrawn from tho fluid.) 

For p =: T, the dynamic equation (appendix B) is 



and therefore 




But 



so that 





i.if 

2 T\ 


( 
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vhere 


T - 





^ - 




Finfally, 


v(c) =■• 


0.332« 




["^ 2 At' 

2 T 


B 


■% 


,=c 


The req.i-iired values of 
table in appendix 3. 


’ © ’ 


and i are obtained from the 


The sinall correction to the slope l(c) is easily calctxlated 
once the mean velocity profile has been obtained (appendix B) . 

Thus 


1(c) 


0.332 n_ _ ^ 

T^ c 


1 ( 0 ) = 0 


The quantity (l - 2?.)v(c) has been calculated as a faction 
of c for various values of T^ at 0, 0.70^ 1.30; 1.50, 2.00 

3.00; and 5.00, and the results of these calculations are given in 
the folloving table. The decisive stabilizing influence of vith- 
draving heat from the fluid at supersonic velocities is illiustrated 
in figure J. 
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1 

Mo = ° 

0.70 

0.1945 

3650 

.80 

,2695 

1080 

.90 

.3485 

402 

1.25 

.5435 

67 

1.50 

.6240 

36 

= I.30; c> 0.231 

0.90 

0.2455 

9230 

1.05 

.4075 

392 

1.20 

.5170 

121 

1.3422 

.5450 

92 

1.50 

.6355 

42 

= 2 . 00 ; c >0.500 

1.63 

0.5074 

671 

1.65 

.5438 

207 

1.70 

.6155 

75 

1.75 

.6749 

40 

1.81 

.7275 

25 

1.85 

.7612 

19 

Mo - ' 

5. 00; c > 0.800 

5.19 

0.8008 

174 

5.20 

.8036 

80 

5.30 

.8262 

23 

5.75 

.9008 

6 

6.0625 

.9350 

3 


^1 

^0 

% 

11 

0 

0 

0.70 

0.1670 

8440 

.80 

.2390 

2110 

.90 

.3265 

613 

1.25 

.5425 

74 

1.50 

.6265 

38 

= I.50; c> 0.333 

1.30 

0.3450 

2770 

1.35 

.4585 

275 

l.ijO 

.5505 

99 

1.4556 

.6276 

49 

1.60 

.7732 

16 

Mp 3. 00; c > 0.667 

2.48 

0.6730 

186 

2.52 

.7058 

59 

2.62 

.7655 

24 

2.72 

.8105 

14 

2.77 

.8295 

10 

2.8225 

.8500 

9 
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APPENDIX D 


V 


SEEWIOR OF 



FROM EQUATIONS OF mm Monow 


In order to determine the effect of free -stream Ifech nimiter, 
thermal conditions at the solid surface^ or f roe-stream pressure 
gradient on laminar stability, it is necessary to know the relation 
between these physical parameters and the distribution of the 


dw ^ 

quantity p — across the boundary layer. The value of — 

a.y 


dw\ 




■J 


at the solid surface is obtained directly from the .i^^namic eciiation 

, ■ 

(equations (6.3) and ( 7 . 2 )). The value of -- 





at the 


sxirfac©^ Tiiich is also useful in the discvinsion of laminar sta- 
bility, is obtained from the dynamic and energy eauations as 
follovrs: 



dv 

3 y 


^ - 


Jl 

.dy^Vl/ 


vj"' 2 vj”T 2/ ''^I’h” ^■^'i'(-i')^ 



Differentiatinc the djaiamic equation once yields the result 


V * ’ • 

^1 




T-, 


- V 


m(m - 1 ) 




+ m 


T-,' 


^1 
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dT* 

At the solid surface the rate of change of ten5>erature 

at* 

— ^p* 

and the rate at vhich the work is done hy pressure gradient u* -- — 

"both vanish, and the rate at 'vdiich a fluid element loses heat by 
conduction equals the rate at "which mechanical energy is transformed 
into heat hy viscous dissipation. The energy equation becomes 



or 



Utilizing the expression for w-^'*' and gives 



- 2 (m + 1 ) 




+ cr(i m )(7 - 1 )Mq^ 



where 




m 4 1 



^I'^l’ 


1 

T ^ 

1 
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From this expression for 


- p / 
d" / v' 


L^y' 


.2 V T 


the following con- 


clusions, which are utilized in the stability analysis, are reached 
When 




vanishes, the quantity 

\ X / 

' /h 

still positive, 

T'Jhen the free- stream velocity is un5.fona. 


d^ /v'^ 


.dj' 


,^Vt 


13 


-Jl 


i(?) 


d 

§7' 


= 0(1 + m)(r - 1)M ^ [ z i ) . + 0(1 ^j2 


M 


.\2 


•P 


3 

V 


\r f • 
''1 } 


that is. 


'—(A 

iy^ /. 


is aliiays positive. 


Idien the surface is insulated. 


d^ / w’ 


•iy" 


u- 


= a(l + m)(7 - 1 )Mq^ 


h'? 




and 


d^ 

_dy^ / 


is always positive, regardless of the pressure 


gradient. 
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APPENDIX E 


a'yLCULATION OP CPITICAL mCH KEMBER POE STABILIZATION 


OF LAimLAE BOUNDAEY lAYM 


For thermal eq.Toili'brluia the rate of heat conduction from the 
gas to the solid surface oalances the rate at vhich heat is radiated 
from the surface. If the rate at which heat is ’'/ri.thdrawn from the 
flviid reaches or exceeds a certain critical valtie at a given local 
supersonic Mach mmber, the laminar boundary- layer fJ.ow is stable 
at all EejTiolds numbers. (See section 6b.) The purpose of the 
following brief calculation is to determine the eciuilibrium surface 
temperatures at several Mach numbers and coii^sare these temperatures 
with the critical temperatures for laminar stability. (See fig. 8. ) 

Nlien the solid surface is in thermal equiilibrium 


where € is the emissivity, a is the Boltzmann constant, and the 
other symbols have already been defined. (See references lii- and I 5 .) 
Consider the case in idiich the free stream is uniform and the 
temperature is constant along the surface. For a = 1, 



iL 


€0 ( 1 ) 


Jo 
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is employed. (See 



O’ 

^1 


if the approximntion p = T 


appendix B . ) Since k.^- 




•d'^O 



0.332 





■When the integrations in eqijation (l) are carried out^ the fol- 
lowing relation is obtained for the detei’mination of the equilibrium 
surface temperature: 


whore 


\|K (t,'i - 1) . (t^ . ^ 



The equilibrium swface teJig)erature under fi*ee-flif^it Rendi- 
tions is affected principally by the vazdation in density with 

altitude h. The results of calculations carried out for alti- 
tudes of 50^000 and 100,000 feet are given in the follo^ving table: 


h 

(ft) 

M 

0 

^3 - ^1 
° -^equil 

® ^cr 

(fig. 8) 

50 X 10^ 

3.0 

0,370 

0.355 

100 X'lo 3 

2.0 

.220 

.185 
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In these calculetione the folloving data are used; 

= 0.50 

L = 2 tb 

T~* = 400° F ahs. 
c 

"a = 4.80 X 10 "^^ Btu/sec/ft^/(de3 F abs,)*^ 

c = 7,73 3 tu/slng/deg F ahs. 

P 

= 3.02 X 10""^ slugs/ft-sec 
a^- = 980 ft /sec 

= 3.61 X 10'^ slugs/ft^ at 50,000 ft 

= 3.31 X 10"^ slugs/ft^ at 100,000 ft 
K = 3.35 X 10'^ at 50,000 ft 
= 3.66 X 10"3 at 100,000 ft 

Since T - Tn > Tg - “ 3 at 50,000 feet 

® equil 

altitude and for = 2 at 100.000 feet a 3 .titude, the laminar 
hoxuadary layer is con^jletely stable under these ccnditions. 

It shovdd be noted that under vind-tunnel-test conditions 3 ji 
which the model is stationary, these radiation- condxrction effects 
are absent, not only because of reradiation fx-om the -walls of the 
•vri.nd tunnel but also because the suxcface temperatures are low - 
generally of the order of room tomperaturo. 
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TABLE I 

AUXILIARY lUNCTIOIIS FOR CALCULATIRG THE STABILITY OF THE LAMUIAR 
BOUNDARY LAYER FOR INSULATED SURFACE 


c 

X 

V 

L 

Hi 

H2 

”3 

N 

3 




Mo 

. 0 




0.0372 

0.0000 

O.OOOl* 

0.0102 

0.5220 

0.2889 

0.0689 

0,2999 

.07UU 

.0001 

.0029 

.0285 

. 14-748 

.2740 

.0604 

.3064 

.1115 

.0003 

.0099 

,0561 

.4303 

.2590 

.0530 

.3124 

.1486 

.0006 

.0235 

,091*0 

.3887 

.2433 

.0460 

.3161 

.1857 

.0012 

.01*62 

.1430 

.3499 

.2278 

.0403 

.3211 

.2226 

.0021 

.0802 

.2040 

.3139 

.2120 

.0350 

.3230 

. 259 ‘^ 

.0033 

. 1281 * 

.2782 

.2808 

.1958 

.0301 

.3217 

. 29 ^ 

.0050 

.1937 

.3670 

.2505 

.1797 

.0256 

.3174 

.3323 

.0071 

.2791* 

.4721 

.2232 

.1639 

.0217 

.3084 

.3682 

.0098 

.3896 

.5960 

.1987 

.11*87 

.0180 

.2935 

.*^037 

.0131 

.5286 

.7418 

.1770 

.1350 

.0139 

.2708 

.im ^3 

.011*2 

.5767 

.7904 

.1711 

.1312 

.0125 

.2618 




Mo - 

0.50 




0.0362 

-0.0000 

- 0.0004 

- 0.0148 

0.5122 

0.2223 

0.0443 

0.1927 

.0723 

-.0000 

-.0001 

-,0234 

.4671 

.2127 

.0401 

.2086 

.1085 

,0001 

.0029 

-.0244 

.4246 

.2019 

.0356 

.2193 

. 14 U 6 

.0003 

.0107 

-.0169 

.3847 

,1904 

.0316 

.2280 

.1806 

.0007 

.0254 

-.0003 

.3474 

,1789 

.0282 

.2366 

.2166 

,0011* 

.0492 

,0260 

.3127 

.1662 

.0249 

.2420 

.2525 

.0023 

. 081*6 

.0627 

,2807 

,1530 

.0217 

,2425 

,2882 

.0036 

,1342 

,1103 

.2513 

.1390 

.0188 

.2406 

.3237 

.0051* 

.2010 

.1695 

,2246 

.1247 

.0158 

.2333 

.3588 

.0076 

.2882 

.2412 

.2005 

.1104 

.0128 

.2179 

.3936 

.0103 

.4000 

.3261 

.1790 

.0963 

.0094 

.1914 

.lt280 

.0137 

.5407 

,4247 

,1602 

.0828 

.0055 

.1444 

.1^306 

.011*0 

.5526 

,4327 

.1589 

.0816 

.0051 

.1397 

. 1*362 

.0146 

.5794 

,4501 

.1560 

.0792 

.0038 

.1262 




Mo- 

. 0.70 




0.0353 

- 0.0000 

-0,0009 

-0.0321 

0.5031 

0,1839 

0.0321 

0.1484 

,0705 

-.0000 

-.0024 

-.0590 

.4599 

.1786 

.0300 

.1652 

.1058 

-.0000 

-.0025 

-.0791 

.4191 

.1721 

,0279 

.1819 

. 11*10 

.0001 

,0006 

-.0914 

.3808 

,1652 

.0257 

.1981 

.1762 

,0004 

.0090 

-.0951 

,341*8 

.1569 

.0233 

.2128 

. 2111 * 

.0008 

.0248 

-.0896 

.3113 

.1478 

.0209 

.2259 

, 21 * 61 * 

.0015 

,0501 

-.0741 

.2802 

.1379 

,0187 

.2358 

.2813 

.0026 

,0872 

-.0478 

.2516 

.1272 

.0165 

,2436 

,3161 

.0039 

.1389 

-.0098 

.2255 

.1157 

,0142 

.2466 

.3505 

.0058 

.2082 

.0412 

,2018 

.1042 

.ou 8 

.2417 

. 381*7 

,0081 

.2985 

,1067 

,1806 

.0925 

.0085 

,2272 

.1+185 

.0109 

.4137 

,1886 

.1619 

,0813 

.0052 

.1987 

. 1+352 

.0126 

.4821 

.2363 

.1534 

,0760 

.0030 

,1787 

.1*1*52 

.0137 

.5270 

.2674 

.i 486 

.0733 

.0016 

.1618 

.4559 

,0149 

.5790 

.3027 

.1436 

.0709 

-.0002 

.1575 
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TABLE I - Concluded 

AUXILIARY FUNCTIONS FOR CALCULATINO THE STABILITY OF THE LAMINAR 
BOUNDARY LAYER FOR INSULATED SURFACE - Concluded 


c 

X 

V 

L 

III 

=2 

M3 

N3 




Mo- 

0.90 




0 . 033 ^ 

.0667 

.1001 

.1335 

.1669 

.2002 

.2335 

.2666 

.2997 

.3326 

.3652 

.3976 

.4296 

.4612 

.4636 

.4812 

0.0000 

-.0001 

-.0002 

-.0002 

-.0001 

.0001 

.0006 

.0012 

.0022 

. 003 ^+ 

.0051 

.0072 

.0098 

.0130 

.0132 

.0153 

-0.0015 

-.0047 

-.0082 

-.0102 

-.0090 

-.0029 

.0098 

.0312 

.0634 

.1086 

.1697 

.2496 

.3518 

.4805 

.>♦913 

.5788 

-0.0503 

-.0972 

-.1389 

-.1746 

..2034 

-.2250 

-.2387 

-.2441 

-.2407 

-.2281 

-.2063 

-.1730 

-,1302 

-.0784 

-.0744 

-.0421 

0.4816 

.4422 

.4049 

.3696 

.3365 

.3054 

.2765 

.2497 

.2251 

.2026 

.1823 

.1641 

.1480 

.13>^0 

.1330 

.1261 

0.1303 

.1298 

.1281 

.1253 

.1213 

.U63 

.U03 

.1030 

.0947 

.0855 

.0759 

.0656 

.0560 

.0464 

.0463 

.0418 

0.0180 

.0185 

.0185 

.0182 

.0175 

.0166 

.0157 

.0143 

.0128 

.ouo 

.0090 

.0060 

.0021 

-.0036 

-.0040 

-.0076 

0.0908 

.1133 

.1366 

. 159 ‘^ 

.1825 

.2055 

.2252 

.2439 

.2597 

.2703 

.2674 

.2515 

.2185 

.1431 

.1373 

.1004 




Mo- 

1.10 




0.0990 

.1320 

.1650 

.19^ 

.2309 

.2638 

.2965 

.3292 

.3616 

.3938 

.4246 

.4572 

.4836 

.5104 

-0.0003 

-.0004 

-.0005 

-.0004 

-.0002 

.0002 

.0009 

.0013 

.0031 

.0049 

.0097 

.0098 

.0126 

.0160 

- 0.0140 

-.0206 

-.0255 

-.0272 

-.0232 

-.0125 

.0072 

.0382 

.0829 

.1442 

.2247 

• 3300 
.4407 

.5789 

-0.2037 

-.2630 

-.3166 

-.3640 

-.4049 

-.4396 

-.4680 

-.4906 

-.5086 

-.5239 

-.5516 

-.5675 

-.6112 

-.6875 

0.4026 

.■^682 

.3358 

.3054 

.2770 

.2506 

.2263 

.2040 

.1837 

.1655 

.1498 

1350 

.1245 

.1151 

0.0673 

.0686 

.0683 

.0667 

.0632 

.0581 

.0516 

.0431 

.0333 

.0218 

.0081 

-.0060 

-.0203 

-.0360 

0.0012 

.0038 

.0051 

.0058 

.0064 

.0062 

.0058 

.0047 

.0031 

.0005 

-.0032 

-.0087 

-.0157 

-.0230 

0.0806 

.1068 

.1319 

.1598 

.1864 

.2101 

.2293 

.2416 
.2454 
.2310 
.1834 
0764 
- 0737 
-.2366 




Mo = 

1.30 



0.2541 

.2858 

.3173 

.3488 

.3800 

.4111 

.4418 

.4721 

.5020 

.5072 

.5416 

-0.0008 

-.0005 

.0001 

.0009 

.0021 

.0037 

.0057 

.0083 

. 011 ^ 

.0120 

.0167 

-0.0561 

-.0505 

-.0364 

-. 0 U 7 

.0258 

.0790 

.1508 

.2449 

.3652 

.3893 

.5777 

-0.5982 

-.6508 

-.6987 

-.7430 

-.7856 

-.8300 

-.8834 

-.9608 

-1.0977 

- 1 . 133 >+ 

-1.3074 

0.2487 

.2255 

.2041 

.1845 

.1667 

.1507 

.1366 

.1242 

.1136 

.1119 

.1020 

0.0244 

.0233 

.0183 

.0109 

.0019 

-.0099 

-.0236 

-.0404 

-.0628 

-.0671 

-.0834 

0.0003 

.0016 

.0014 

.0003 

-.0016 

-.0048 

-.0090 

-.0169 

-.0294 

-.0324 

-.•0549 

0.2200 

.2440 

.2644 

.2742 

.2700 

.2285 

.1184 

-.0818 

- . 4943 
-.5971 
-1.5080 
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TABIZ II 

AUXIUAPY FDNCTIOHS TOR CAICTIATIIIIG THE STABIUTI OF THE 
lAMINAB BODHDAHY lAYER FOB NOHIHSUIATED SURFACE 


c 

X 

V 

L 

Hi 

®2 


"3 




- 0.10] 

1 - 0.70 




0.0262 

0.0056 

0.0825 

0.0635 

0.6102 

0.3272 

0,0524 

0.2748 

.0521 

.0112 

.1645 

.0949 

.5725 

.3157 

.0502 

.2920 

.0777 

.0166 

,2k66 

.1184 

.5367 

.3045 

.0481 

.3081 

.1030 

.0220 

.3297 

,1400 

.5026 

.2936 

.0458 

.3233 

.1281 

. 027 *^ 

.U1U6 

.1632 

.4703 

.2828 

.0433 

.3380 

.1529 

.0327 

.5023 

,1904 

.4396 

,2724 

,0412 

.3519 

.1701 

.0365 

.5661 

,2130 

.4191 

.2651 

.0396 

.3610 

.1726 

.0370 

.5754 

.2163 

.4162 

.2642 

.0394 

.3623 




. 0 . 70 ; 

- 0.80 



0.0237 

0,0033 

0.0486 

0.0279 

0.5954 

0.28U 

0.0493 

0.1369 

.0U72 

.0066 

.0965 

.0374 

.5620 

.2737 

.0475 

.1504 

.0705 

.0099 

.1443 

.0430 

.5300 

.2663 

.0457 

.1635 

.0937 

.0132 

.1925 

,0482 

,4994 

.2590 

.0437 

.1763 

.1168 

.0164 

.2417 

.0550 

.4701 

.2514 

.0417 

.1882 

.1397 

.0197 

.2926 

.06^9 

.4420 

.2439 

.0397 

.2001 

.1625 

.0230 

.3457 

.0789 

,4152 

.2363 

.0378 

.2110 

.1851 

.0263 

.4017 

.0982 

.3897 

.2287 

.0359 

.2213 

.2075 

.0297 

,4614 

.1236 

.3654 

.2210 

.0339 

.2311 

.2298 

.0331 

.5253 

.1562 

.3424 

.2133 

.0321 

.2400 

. 2 »t 09 

. 03 V 9 

.5592 

.1754 

.3313 

.2094 

.0310 

.2443 

. 2»*75 

.0359 

,5801 

.1877 

.3248 

.2071 

.0303 

.2465 




- 0 . 70 ; 

T^^ - 0.90 




O.OU33 

0.0036 

0.0517 

0.0051 

0.5506 

0.2410 

0.0435 

0.1426 

.0863 

.0072 

.1028 

-.0047 

.4939 

.2304 

.0404 

.1638 

.1291 

.0108 

.1568 

-.0111 

.4414 

.2191 

.0370 

.1846 

.nih 

.0145 

.2173 

-.0079 

.3930 

.2074 

.0337 

.2032 

.2135 

.0185 

.2885 

.0096 

.3485 

,1951 

.0304 

.2203 

.2551 

.0227 

.3746 

,0462 

.3080 

.1825 

.0272 

.2339 

.2963 

.0274 

.4805 

.1073 

.2715 

.1696 

.0240 

.2462 

.3166 

.0299 

. 5 *t 26 

.1489 

.2547 

.1637 

.0224 

.2517 

.3268 

.0312 

.5762 

.1776 

.2466 

.1606 

.0217 

.2541 




Mj, - 0 . 70 ; 

Ti - 1.25 




0 . 03 ^ 

-0.0016 

-0.0237 

-0.0476 

0.5100 

0.1750 

0.0324 

0.1462 

.0692 

-.0032 

-.0476 

-.0797 

.4678 

.1710 

.0310 

.1634 

.loUo 

-.0048 

-.0698 

-.1013 

.4276 

.1661 

.0292 

.1794 

.1389 

-.0062 

-.0886 

-.1132 

.3896 

.1600 

.0272 

.1956 

.1738 

-.0076 

-.1021 

-.1155 

.3538 

.1529 

.0251 

.2108 

.2088 

-.0087 

-.1085 

-.1081 

.3202 

.1448 

.0228 

.2238 

. 2 U 39 

-.0095 

-.1057 

-.0912 

.2888 

.1354 

.0208 

.2342 

.2789 

-.0101 

-.0917 

-.0645 

.2597 

.1249 

.0185 

.2402 

.3138 

-.0103 

-.0641 

-.0281 

.2330 

.1133 

.0161 

.2409 

. 3 >t 85 

-.0100 

-.0203 

,0179 

.2086 

.1008 

.0139 

.2297 

.3831 

-.0092 

.0427 

.0734 

.1865 

.0870 

.0113 

.2069 


-.0079 

.1286 

.1373 

.1668 

.0728 

.0083 

.1616 

.U512 

-.0059 

.2414 

.2071 

.1495 

.0582 

.0042 

.0816 

.U 846 

-.0031 

.3859 

.2770 

.1345 

.0427 

-.0012 

-.0601 

.5092 

-.0006 

.5184 

.3212 

.1248 

.0314 

-.0067 

-.2262 

.5190 

.0006 

.5779 

.3349 

.1212 

.0269 

-.0091 

-.3028 
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TABUE m 

PHASE VELOCITY, WAVE HTJ^BERjAND REYNOLDS NIMBER FOR NEOTRAL SUBSONIC 
DISTURBANCE (STABHITy LDGTS) FOR INSaLATED SURFACE 




0.0321 

25,500,000 

.0685 

1,500,000 

.1103 

278,000 

.1585 

83,000 

.2146 

32,600 

.2808 

14, 800 

.3590 

7,700 

.4535 

4,420 

.5707 

2,760 

.7243 

1,850 

.9589 

1,360 

1.0770 

1,280 

1.2730 

1,530 

1.2940 

1,880 

1.1980 

3,530 

1.0400 

6,710 

.8728 

13,300 

.7177 

27,500 


0.0372 

.0744 

.1115 

.1486 

.1857 

.2226 

.2594 

.2960 

.3323 

.3682 

.4037 

.4143 

.4143 

.4037 

.3682 

.3323 

.2960 

.2594 


0.0362 

.0723 

.1085 

.1446 
.1806 
.2166 
.2525 
.2882 
.3237 
.3588 
.3936 
.4280 
.4306 
.4362 
.4362 
. 4306 
.4280 
3936 
.3588 
.3237 
.2882 


0.0038 

.0082 

.0131 

.0189 

.0255 

.0334 

.0427 

.0540 

.0679 

.0862 

.1142 

.1282 

.1515 

.1540 

.1424 

.1238 

.1039 

.0854 


0.0251 

35,600,000 

.0538 

2,130,000 

.0868 

392,000 

.1250 

116,000 

.1695 

44,500 

.2216 

20,200 

.2829 

10,lt00 

.3556 

5,850 

.4442 

3,570 

.5549 

2,330 

.6993 

1,620 

..9301 

1,230 

.9558 

1,220 

1.0140 

1,190 

1.1880 

1,410 

1.2150 

1,580 

1.2150 

1,660 

1.1240 

3,080 

.9788 

5,670 

.8272 

10,800 

.6869 

21,100 


Mo - 0.50 


3,030,000 

178,000 

33,100 

9.880 

3.880 
1,760 

917 

526 

329 

220 

162 

153 

182 

223 

421 

799 

1,580 

3,270 


0.0029 

4,270,000 

.0063 

248,000 

.0101 

45,700 

0146 

13,500 

.0198 

5,190 

.0258 

2,360 

.0330 

1,210 

.0414 

682 

.0518 

4l6 

.0647 

272 

.0815 

189 

.1084 

144 

.1114 

142 

.1182 

139 

.1384 

164 

.I4l6 

184 

.1416 

194 

.1310 


.1141 

661 

.0964 

1,260 

.0800 

2,460 


0.0353 

.0705 

.1058 

.1410 

.1762 

.2114 

.2464 

.2813 

.3161 

.3505 

.3847 

.4185 

4352 

.4452 

.4559 

.4559 

.4452 

.4352 

.4185 

. 3847 

.3505 

.3161 


0.70 


0.0191 

53,400,000 

.0415 

3,060,000 

.0677 

555,000 

.09^ 

161,000 

.1344 

61,100 

.1766 

27,300 

.2268 

i3,800 

.2857 

7,630 

.3570 

4,550 

.4433 

2,900 

.5‘=>15 

1,960 

.6951 

1,420 

.7917 

1,230 

.8655 

1,160 

.9704 

1,110 

l.i230 

1,330 

1.1420 

1,650 

1.1230 

1,930 

1.0720 

2,670 

.9"^81 

.7965 

4,810 

8,880 

.6659 

16,700 


0.0022 

6.100,000 

.004? 

349,000 

.0077 

63,400 

.0112 

18,400 

.0154 

6,980 

.0202 

3,120 

.0259 

1,580 

.0326 

872 

.0408 

520 

.0506 

331 

.0630 

224 

.0794 

162 

.0904 

141 

.0939 

132 

1108 

127 

.1283 

152 

1304 

189 

.1283 

227 

.1225 

305 

.1072 

550 

.0910 

1,010 

.0761 

1,910 
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TABLE III - Concluded 

PH.^SE VELOCITY, WAVE MMBER, AND REYNOLDS NUMBER FOR NEUTRAL SUBSONIC DISTURBMCE 
(STABILnY LIMITS) FOR INSULATED SURFACE - Concluded 


c 

a 

R 




Mq - 0.90 

0.0334 

0,0107 

111,000,000 

0.0012 

12,600,000 

.0667 

.0248 

5,960,000 

.0028 

679,000 

.1001 

.0421 

1,030,000 

.0048 

117,000 

.1335 

.0632 

290,000 

.0072 

33,000 

.1669 

.0885 

106,000 

.0101 

12,100 

.2002 

.1186 

46,000 

.0135 

5,240 

.2335 

.1540 

22,600 

.0175 

2,570 

.2666 

.1961 

12,100 

.0223 


.2997 

.2459 

7,020 

.0280 

800 

.3326 

.3053 

4,320 

.0348 

492 

.3652 

.3777 

2,820 

.0430 

321 

.3976 

.4638 

1,950 

.0529 

222 

.4296 

.5733 

1,410 

.0653 

161 

.4612 

.7260 

1,090 

.0827 

125 

.4636 

.7396 

1,080 

.0843 

123 

.4812 

.8810 

1,010 

.1004 

115 

. 4812 

1.0130 

1,230 

.1154 

140 

. 4636 

1.0120 

1,740 

.1153 

199 

.4612 

1.0070 

1,820 

.1148 

207 

.4296 

.9027 

3,180 

.1029 

363 

.3976 

.7823 

5,590 

.0892 

637 

.3652 

.6642 

9,940 

.0757 

1,130 

.3326 

.5523 

18,500 

.0629 

2,100 


0.0990 

.1320 

.1650 

.1980 

.2309 

.2638 

.2Q6^ 

.3292 

.3616 

.3938 

.k2h6 

A372 

.1^36 

.^lOh 

.^lOk 

.4836 

. 45.72 

.4246 

.3938 


Mo . 1.10 


0.0086 

5,730,000 

.0268 

769,000 

.0468 

224,000 

.0707 

85,000 

.0991 

38,300 

.1329 

19,300 

.1727 

10,600 

.2200 

6,260 

.2755 

3,920 

.3417 

2,610 

.4159 

1,850 

.5193 

1,350 

.6268 

1,100 

.8010 

991 

.9165 

1,220 

.8927 

2,060 

.8023 

3,320 

.6785 

5,930 

.5766 

10,400 


Mo 


1.30 


0.0009 

618,000 

.0029 

82,900 

.0050 

24,100 

.0076 

9,160 

.0107 

4,130 

.0143 

2,080 

.0186 

1,140 

.0237 

675 

.0297 

423 

.0368 

281 

.0448 

199 

.0560 

146 

.0676 

119 

.0864 

107 

.0988 

131 

.0962 

223 

.0865 

3-58 

.0732 

639 

.0622 

1,120 


0.2541 

.2858 

.3173 

.3488 

.3800 

4111 

.4418 

.4721 

.5020 

.5072 

.5416 

.5416 

.5072 

.5020 

4721 

.4418 


0.0451 

.0818 

.1202 

.1636 

.2132 

.2707 

.3377 

.4166 

.5123 

.5316 

.7582 

.8931 

.7781 

.7592 

.6458 

.5401 


63.800 

24.800 

12, 300 
6,990 
4,280 
2,800 
1,930 

1,420 

1,110 

1,070 

886 

1,080 

2,310 

2,550 

4,500 

7,980 


0.0047 

6,630 

.0085 

2,570 

.0125 

1,280 

.0170 

726 

.0222 

445 

.0281 

291 

.0351 

201 

.0433 

147 

.0532 

115 

.0552 

111 

.0788 

92 

.0928 

112 

.0809 

241 

.0789 

265 

.0671 

468 

.0561 

829 
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TABLE IV 

% 

PHASE VELOCITY, WAVE NIMBER, AND REYNOLDS NIMBER FOR NEUTRAL 
SUBSONIC DISTURBANCE (STABILITY LIMITS) FOR NONINSULATED SURFACES 


c 

» a 

R 



> 0 . 70 ; Tx » 0.70 

0.0262 

.0521 

.0777 

.1030 

.1281 

. 1'?29 

.1701 

.1726 

.1726 

.1701 

.1529 

,1281 

.1030 

.0777 

.0521 

0,0339 

. 073 ^^ 

,1188 

.1708 

,2308 

.3030 

.3670 

.3777 

,1+986 

. 1+977 

,1+732 

. 1+175 

,3^60 

.2620 

.1713 

82,1+00,000 

5.360.000 

1.110.000 

371.000 

161.000 
83,1+00 
57,200 

5 l+,l +00 

69,000 

73,900 

121.000 

270.000 

711.000 

2.500.000 
ll+, 600,000 

0 004 l 
.0088 

. 011+3 

.0205 

.0277 

.0361+ 

, 0 l+ 4 l 

.0454 

.0599 

.0598 

.0568 

.0502 

.0416 

.0315 

.0206 

9.900.000 

644.000 

133.000 
44,600 
19,300 
10,000 

6.870 
6,540 
8,280 

8.870 

14,500 

32.400 

85.400 

300.000 

1.750.000 

Mq - O.7O; T^ - 0.80 

0.0237 

. 01^72 

.0705 

.0937 

.1168 

.1397 

.1625 

.1891 

.2075 

.2298 

. 2 i +09 

. 21+75 

, 21+75 

,21+09 

,2298 

.2075 

.1851 

, 162'5 

.1397 

.1168 

.0937 

.0705 

0.0237 

.0501+ 

.08ol+ 

.1138 

.1509 

.1923 

.2382 

.2908 

.3510 

. 1+237 

,1+668 

,1+962 

.6308 

,6233 

.6056 

,5609 

.5062 

.1+1+69 

.3827 

.3161+ 

,21+89 

.1822 

157,000,000 

9.910.000 

1.970.000 

633.000 

263.000 

129.000 

70.900 
1+2,600 

27.500 
18,800 

15.900 

ll +,500 

18.500 
21,1+00 
27,200 
l+l+,900 
77 , 1+00 

11+1,000 

280.000 

630.000 

1.690.000 
5,8^,000 

0.0028 

.0059 

.0094 

.0133 

.0176 

.0224 

.0278 

.0339 

.0409 

.0494 

.0544 

.0578 

.0735 

.0726 

.0706 

0654 

.0590 

.0520 

.0446 

.0369 

.0290 

.0212 

18,300,000 

1,150,000 

230.000 
73,700 

30.600 
15,000 

8,260 

4,960 

3,200 

2,190 

1,860 

1,690 

2,160 

2,500 

3,170 

5,230 

9,010 

16. 400 

32.600 

73.400 

197.000 

686.000 
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TABLE 17 - Concluded 

PHASE VELOCITY, VrAVE NUMBER, AND EEYNOLDS NUMBER FOE 
NEUTRAL SUBSONIC DISTURBANCE (STABILITY LIMITS) FOE 
NONINSULATED SURFACES - Concluded 


c 

a 

R 





Mq . 0.70; Tj^ - 

0.90 


0,0433 

.0863 

.1291 

.1714 

.2135 

.2551 

.2963 

.3166 

.3268 

.3268 

.3166 

.2963 

.2551 

.2135 

.1714 

.1291 

,0863 

0.0368 

.0815 

.1353 

.1996 

.2775 

.3728 

.4980 

.5814 

.6347 

.7817 

.7701 

.7307 

.6275 

.5133 

.3972 

.2858 

.1793 

17,100,000 

1 . 040.000 
200,000 

62.500 

25.500 
12,400 

6,970 

5,520 

4,990 

6,500 

7,920 

11,600 

25,200 

60,300 

170.000 

617.000 

3.740.000 

0.0042 

.0092 

.0153 

.0226 

.0314 

.0422 

.0563 

.0658 

.0718 

.0884 

.0871 

.0827 

.0710 

.0581 

.0449 

.0323 

.0203 

1,930,000 

118.000 
22,700 

7.070 

2,880 

1,410 

789 

624 

565 

735 

895 

1,310 

2,850 

6,820 

19,200 

69,800 

423.000 



Mq =1 0 . 70 ; T^ > 

1.25 


0.0346 
. 0692 
.1040 

.1389 

.1738 

.2088 

.2439 

.2789 

.3138 

.3485 

.3831 

. 4174 

.4512 

.4846 

.5092 

.■5190 

.5190 

.5092 

. 4846 
.4512 
.4174 

0.0160 

.0346 
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Fig. 1 
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Figure !•- Boundary-layer velocity profiles for insulated surface* 
Since 6 is taken equal to unity, the temperature profile is 


given hy T = 


. [ Tj - (l 4 
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Fig. 2 



Figure 2.- Boundary-layer velocity profiles for nonlnsulated 
surface. = O.70. Is the ratio of surface temperature 

(deg abs.) to free-stream temperature (deg abs.). Inflection 
more pronounced and farther out Into fluid for =* 1«25 

than for Insulated surface (T]^ = 1.10). No Inflection for 

0.70, 0 .' 60 , 0.90. 
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(a) Inaulatod surface. 

Figure 3.- Diatrlbutlon of across boundary layer, as 

expressed by the function ¥( 0 ) which appears directly In 
stability calculation. The value of c ■ c© at which 
(1 - 2X)v ■ 0,580 ( \ is small) Is a measure of the 
stability of a given laminar boundary^-layer flow, 

R. « 6 . 

t7 


'cr 


min ^ Cq^ /l - ^f(l - Oq)^ 


(,.)3dyiT4 
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Fig. 3b 


O 




I 


o 



c 

(1)) Noninsulated surface. Mq s 0«70. 


Figure 3«- Concluded^ 
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(a) 

rigur« 4.- 


Insulated surface* 
Ifore nuaber ocq 


• 0 (Blasius flow)* 

against Reynolds nuaber R© for neutral stability 
of laalnar boundary- layer flow. 
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Cb) Insulated surface. Mq = 0.50. 


Figure 4.- Continued. 


NACA TN No. 1360 Fig. 4b 




(c) Insulated surface. Hq = 0.70* 
Figure 4«* Continued. 


Fig. 4c NACA TN No. 1360 




(d) Insulated surface* Mq e 0*90* 
Figure 4*- Continued# 
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Fig. 4e NACA TN No. 1360 
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(f) Insulated •urfAoa* 
Figure Contixmed. 


= 1 ^ 30 . 


NACA TN No. 1360 Fig. 4f 




Fig. 4g 
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(g) Monlnfrulatad stirface. Mq * 0.70j T^l * (~) >0 

^07/j^q 

Figure 4*- Continued# 
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(1) Voalnvulated mirfsioe. s 0«70| 
Figure 4«- Continued* 


0,90j /2I) > 0. 

'®T/y =0 


Fig. 4 i NACA TN No. 1360 
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(k) Insulated surface. Note that a-> Og 0 as R-> » for insulated surface. (Mq ^ 0.) 
Figure 4 . - Continued. 
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(/) Noninsulated surface. Mq = 0.70. Note that 0 as R-> « when = 1.25. 

Figure 4.- Concluded. 
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Fig. 5 
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Figure 5*- Variation of minimum critical Reynolds ninnber 
with Maoh number for laminar boundary-layer flow. 

/ u ~ ^0 \ 

Insulated surface. Rq = 1 > (“^4^ — i 
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Fig. 6a 
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figure 6.- Dependence of minimum critical Reynolds number 
on thermal conditions at solid surface. = O.70. 

Is the ratio of surface temperature (deg abs.) to 
free-stream temperature (deg abs.). 


Fig. 6b 
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(b) ^ against oaloulated from Taluee 

*crmln ^ 


of R taken from figure 4. 

^ornin 


Figure 6.- Oonoluded 
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Fig. 7 
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Figure 7,- Stabilizing effect at supersonic Mach numbers of 
withdrawal of heat from fluid through solid surface. At 
each value of M > 1, there Is a critical value of s 


such that for Tt ^ Tt the laminar boundary-layer flow Is 

i Xqj. 


stable at all values of the Reynolds number. (Curves for 

M = 0 and Mrt = O.yO Included for comi>arl8on. ) Rq 
o o ®cr„ 4 « 

estimated from equation CT s 1. 


or 
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Fig. 8 
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Figure S.- Critical temperature ratio 


for 


or 


stability of laminar boundary layer against 
Mach number M^, Is the ratio of stagnation 

temperature (deg abs*) to free-stream teirperature 


L.^ 


for cr = 1 
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Fig. 9 
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